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ABSTRACT 

We study the stellar populations of Type la Supernovae (SNe la) host galaxies using 
Sloan Digital Sky Survey (SDSS)-II spectroscopy. We focus on the relationships of 
^ , SNe la properties with stellar velocity dispersion and the stellar population parameters 

age, metallicity and element abundance ratios. We utilize stellar population models 
of absorption line indices for deriving the stellar population parameters. Furthermore, 
we revisit the correlation between Hubble residual and photometric galaxy stellar 
mass, which has received much attention in the recent literature. We concentrate 
on a sub-sample of 84 SNe la from the SDSS-II Supernova Survey with available 
host galaxy spectroscopy. This sub-sample has been selected based upon the quality 
of the spectroscopy and accuracy of the derived SNe la parameters. In agreement 
with previous findings, we find that SALT2 stretch factor values show the strongest 
dependence on stellar population age in terms of a clear anti-correlation. Hence, SNe la 
peak-luminosity is closely related to the age of the stellar progenitor systems, where 
r> , more luminous SNe la appear in younger stellar populations. We find no statistically 

^ ■ significant trends in the Hubble residual with any of the stellar population parameters 

studied, including age and metallicity contrary to the literature, as well as with stellar 
velocity dispersion. We extend the sample to also include SNe la with available SDSS 
host galaxy photometry only. We find that the method of stellar mass derivation 
is affecting the Hubble residual-mass relationship when lower number statistics are 
used, while this effect is weaker for the extended sample. For this larger sample (247 
objects) the reported Hubble residual-mass relation is strongly dependent on the stellar 
mass range studied and behaves as a step function. In the high mass regime, probed 
by our host spectroscopy sample, the relation between Hubble residual and stellar 
mass is flat. Below a stellar mass of ~ 2 x 10 10 Af Q , i.e. close to the evolutionary 
transition mass of low-redshift galaxies reported in the literature, the trend changes 
dramatically such that lower mass galaxies possess lower luminosity SNe la after light- 
curve corrections. This non-linear behaviour of the Hubble residual-mass relationship 
should be accounted for when using stellar mass as a further parameter for minimising 
the Hubble residuals. 

Key words: supernovae: general cosmology: observations distance scale cosmolog- 
ical parameters large-scale structure of Universe - galaxies: abundances 
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1 INTRODUCTION 

Type la supernovae (SNe la) are useful for constraining 
cosmological parameters. Their large peak luminosities can 
probe vast cosmological distances and connect redshift space 
to luminosity distance. This property led t o the discovery of 
an accelerating expans ion of the universe (|Riess et al.ll 199^ ; 
iPerlmutter et alJll999h . 

To be ideal cosmology indicators, all SNe la explosions 
should have the same peak luminosity, but this is not the 
case. Observed SNe la span a range in peak lumin osity 
accompanied by varying decline rates (|Phillips! Il993l ). i.e. 
the peak luminosity decreases with increasing decline rate. 
This variation can be corrected for to find a standardised 
peak luminosity. Several light-curve fitting tools are avail- 
able djha et all 120071 : iGuv et all 120071 ; IConlev et all |200S| ; 
iKessler et al.ll2009bh . where the shape of the light-curve and 
the colours are corrected to match a standardised peak lumi- 
nosity. The light-curve shape correction is known as stretch- 
factor. 

Correcting the light-curves increases the precision of 
the derived luminosity distances and consequently reduces 
the scatter in the redshift-distance relation, thus increasing 
the precision of the derived cosmological parameters. How- 
ever, even after light-curve corrections this scatter is non- 
negligible. Understanding systematic uncertainties in the de- 
rived SNe la light-curve parameters is therefore key to im- 
proving supernovae cosmology. 

The accepted model for SNe la is thermonuclear 
explosion of a carbon-oxygen whi te dwarf (WD) that 
reaches the Chandrasek h ar lim it (|Whelan fe Ibenl 1 19731 ; 
iHillebrandt fe Niemeveil l200d ). Two different chan- 
nels have been proposed, either the single-degenerate 
(SD) scenario where mass is accreted from an evolved 
main-sequence binary companion, or the double- 
degenerate (DP) c ase o f merging of t wo WDs (e.g 



19861; IBranch et al.l 



Woosley fc Weaverl 

2001; Hoflich, ~Kh"okhlov fc Wheeler! [l995; 



1995; IBranch! 



Grcggio 2005: 



Yungelson fc Livid 120001 )" This raises the possibility that 
different SNe la populations may be present. On the 
other hand, it has been shown that a DD system is likely 
to lead to an accretion -induced collapse rat her than a 
thermonuclear explosion (|Saio fc Nomotdll998h . 

The delay time, i.e., the time between progenitor for- 
mation and explosion, in the DD scenario is determined by 
the life-time of the WD progenitors and the orbit of the 
two binary stars. The delay-time of a SD system partly de- 
pends upon the main-sequence lifetime of the companion 
star. A wide variety of delay times h ave been observation 
ally suggested, rangi ng from <1 G vr dBarris fc Tonrvl 
Aubourg et~ai1l2008ll to > 2 Gyr (|Gal-Yam fc Maod 



2006 



2004 



Strolger et al 1l2004l2005h . Considering several different pro 



genitor systems, theoretical models find the rate of SNe la 
explosi ons (SNR) to peak at delay times below or close to 
1 Gy r (|Yungelson fc Lividl2000l ; lGreggidl2005l ; iRuiter et al.l 
2010). The SNR of most progenitor systems then smoothly 
declines and becomes 10-100 times lower at delay times of 
~10 Gyr. Comparing observed delay times to theoretical 
predictions can constrain possible progenitor systems. 

It is well established that the SNR is higher in 
star-forming late-t ype than in passively evolving early- 
type galaxies (e.g. lOemler fc Tinslevl Il979l ; Ivan den Berghl 



199d ; lMannucci et al.l2005l ; ISullivan et a l. 2006; S mith et al.l 



20121 ). Moreover, several authors found a d ependence on host 
galaxy mass for the de c line rate of SNe la (jKellv et allfeoid ; 
ILampeitl et all l2010bl ; ISullivan et~al"1 l20ld ). Since galaxy 
mass correlates with stellar popula tion parameters and 
properties of the inter-stellar medium (flremonti et al.|[20o3 : 
iGallazzi et al.|[2~005l ; iThomas et al.ll2010l ). more fundamental 
correlations such as with stellar population age, metallicity 
and element abundance ratios may be expected. The tight 
relation between SNe la decline rate and peak brightness in- 
dicates a primary dependence of luminosity on one of these 
parameters, thus holding important information about the 
properties of SNe la progenitor systems. 

The luminosity of SNe la arises from the radioac- 
tive decay of 56 Ni to ^ 6 Co, that t hen decays to 56 Fe 
(|Colgate fc McKeel Il969l : lArnettl Il982f). such that the peak 



brightness depends on the 56 Ni mass. iTimmes" et all (|2003h 
show theoretically that metallicity effects, in the range 1/3- 
3 Zq, can induce a 25% variation in the 56 Ni mass. The- 
oretical models have also predicted that the carbon mass 
fraction, determined by the metallicity and mass of the WD 
progenitor, is respo nsible for SNe la luminosity variations 
(|Umeda et al.|[l999r ). Dependencies of SNe la properties on 
element ratios, such as C/Fe, may thus be expected. The 
ratio of light elements to Fe is also interesting because Fe 
dominates the light-curve of SNe la and could well influence 
some of the SN la properties. 

Observationally, it has been found that SNe la in star- 
forming galaxies show slower decline rates as compared 
to S Ne la in passively evolving galaxies (ISullivan et al 



200d; IHowell et all lipo! iNeill et all 120091 ; ILampeitl eta! 



2010d : ISmith et all 120121 ). Moreover, several authors have 
recently found dependencies on stellar p opulation age 



and / o r metallicity for SNe la decline rate (iHamuv et al 



2000 |; IGallagher et alJl200St IHowell et all 120091 ; INeill et al 
200S ; iGupta et al.ll201ll ). However, eit her fairly smal l sam- 
ples have been used (< 30 objects, IHamuv et all l200d 



IGallagher et all 20081) or metallicity has bee n measured in 
directly ([Howell et alj|2009l ; INeill et al.ll2009l ). 

It is therefore desired to study the full range of stellar 
population parameters for a statistically significant sample 
of SNe la host galaxies, which was the main aim of this 
study. 

Correlations between SNe la peak luminosity and stellar 
population parameters should ideally be eliminated through 
light-curve fitting, due to the tight correlation between peak 
luminosity and decline rate. However, if such correlations re- 
main even after light-curve corrections, the previously men- 
tioned non-negligible scatter in the redshift-distance rela- 
tion could be reduced. These relationships are important to 
identify, especially for large redshift surveys, as the stellar 
population parameters age and metallicity change over cos- 
mic time. Since the scatter in the redshift-distance relation 
is small, high-quality data are required to detect dependen- 
cies on host galaxies. To date galaxy mass and metallicity 



(IGallagher et all 


20081; IHowell et al.ll2009l; INeill et all 


200S; 


Kellv et al.ll201fi; 


Lampeitl et al.ll2010d;ISullivan et all 


201C; 


Gupta et al. 2011; D'Andrea et al. 2011). This led to the 



introduction of stellar mass as a further parameter, besides 
light-curve stretch-factor and colour, to minimise the scatte r 
in the redshift-distance relation (e.g. ILampeitl et al.ll2010"bl ). 
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However, it is desirable to find the fundamental parameter 
of this scatter and, again, to do this the full range of stellar 
population parameters is needed for a statistically signifi- 
cant sample. 

In this paper, we study host galaxies of SNe l a from the 
Sloan Digital Sky Survey (SPSS jYork et alJl200Ch -II Super- 
nova Survey |Frieman et al.l l2008l) for which SDSS-II spec- 
troscopy is available. This allows us to investigate SNe la 
host galaxy stellar populations through Lick absorpt ion line 
indices (e.g. IWorthev et al ] |!994lTrager et al.lll998l ). These 
are defined for 25 prominent absorption features in the op- 
tical and are useful for breaking the age-metallicity degen- 
eracy, deriving element ratios an d they are insensi t ive to 
dust reddening |Mac Arthur! 120051 ). iGallagher et~ai] |2008) 
is to date the only study in the literature that relies on 
Lick indices for studying SNe la h ost galaxies, based on 
a sample of 29 early-type galaxies (|Gallagher et al.l 120081 . 
study a sample of five SN la host galaxies through Lick 
indices). The SDSS-II Supernova Survey in combination 
with SDSS-II spectroscopy allow us to significantly im- 
pr ove upon this sample siz e. We use the method presented 
in Ijohansson et al.l (|2012l ). which is based on up-to-date 
stellar population models of ab s orptio n line indices from 
iThomas. Maraston fc Johansson] (|201ll . TMJ), to derive age 
and metallicity together with a range of element abundance 
ratios including C/Fe as well as O/Fe, Mg/Fe, N/Fe, Ca/Fe 
and Ti/Fe. This allows us to explore the full range of stellar 
population parameters to search for the fundamental param- 
eter of the variations of SNe la properties. 

In addition, we determine stellar masses from SDSS- 
II photometry using different methodologies of photometric 
SED-fitting to study systematic uncertainties in the pro- 
posed dependency of the scatter in the redshift-distance re- 
lation on stellar mass. 

The paper is organised as follows. The data sample used 
is described in Section [2] along with the description of our 
derivation of spectroscpic stellar population parameters and 
stellar masses. The relationships between SN la properties 
and host galaxy parameters are presented in Section [3] We 
discuss the results in Section 2] and concluding remarks are 
given in Section [5] 



2 DATA SAMPLE 

The sample used in this study is drawn from the SDSS- 
II Supern ova Survey jFrieman et alJl200Sh. T he SDSS 2.5 m 
telescope (|York et alj200d : lGunn et alj2006t ). located at the 



120091 ; iLampeitl et~aH |2010al ) as well as photometric and 
spectroscopic studies of SNe la host galaxy properties 



spectroscopic studies 01 alMe la nost galaxy properties 
(ILampeitl et all l2010bl; iGupta et alj 1201 if : iD'Andrea et all 
l201ltlsmith et alj|20ia>. 



Apache Point Observ atory (APO), is e quipped with a multi- 
object s pectrograph (ISmee et al.ll2012T i and wide field CCD 
camera (G unn et al.lll9980 . During the eight year period of 
2000-2008, the SDSS-I and SDS S-II surveys obtained deep 
images in the SDSS ugriz filters (|Fukugita et al.lll996l ) and 
spectroscopy for more than 930,000 galaxies. 

The SDSS-II SN survey was performed over a three year 
period (2005-2007) to repeatedly image transient objects in 
the SDSS "Stripe 82" region. Thousands of potential SN 
candidates were observed. Out of these, 890 were identified 
as SN la candidates including 551 with spectroscopic confir- 
mation dSako et al.|[20oi : iHoltzman et alj|200g| ; ISako et all 
I in prepj) . The SDSS -II SN sample has b e en used for cos- 
mological analyses |Kessler et al.l l2009al : ISollerman et al.l 



Following Lampeitl et al.l (|2010bl ), we include photo- 
metrically classified SNe la, to ensure a more co mplete sarn- 
ie. B ased on the Bayesian light-curve fitting of lSako et al.l 
20081) . these supernovae have a light-curve consistent with 
being a Type la. The contamination of non-SNe la objects 
using the photome tric classification scheme is only ~3% 
jPildav et al.ll2010l ). 

In this work we focus on a sub-sample of SNe la with 
available SDSS-II host galaxy spectroscopy (from now on 
referred to as the host spectroscopy sample). The SN la co- 
ordinates and redshifts have been cross-matched with the 
corresponding parameters in the SDSS-II DR7 catalogue to 
identify the nearest object with a spectrum within a 0.25" 
radius, resulting in a host spectroscopy sample of 292 ob- 
jects. A similar cross-match has been made to ac quire co-add 
phot ometry from the SDSS "Stripe 82" region (|Annis et al.l 
l201lf ). The host galaxy coordinates from both runs have 
been further cross-matched to ensure that the spectroscopy 
and photometry belong to the same object. Out of the 292 
supernovae in the host spectroscopy sample, 138 (47%) are 
spectroscopically confirmed SNe la. For redshifts below 0.21, 
the fraction of spectroscopically confirmed SNe la increases 
to 60%. This is also the redshift range covered by the final 
sample selection (see Section r2.3[) used in the main analysis. 

In addition to the host spectroscopy sample we include 
in Section 13.4.21 SNe la host galaxies with available SDSS 
photometry only, to extend the stellar mass analysis. This 
extended ph otometric sample adds a nother 102 objects from 
the study of ILampeitl et all (|2010bl ) and is described in de- 
tail in Section f3. 4. 21 

The derived SNe la properties, host galaxy parameters 
and final selection of the host spectroscopy sample are de- 
scribed in the following sections. 



2.1 SNe la properties 

The lig ht- curve fitting technique follows ILampeitl et al.l 
(|2010bl ) to which the reader is referred for more details. A 
summary of the main features is presented here. 

A number of light-curve fitting methods exist in the lit- 
erature and several authors have confirmed tha t the results 
are r o bust against different analysis algorithms (jlCellv et al.l 
l20ld: ILampeitl et all l2010bl : ISullivan et ail |2010| ). SALT2 
( Guv et al.1 120071 ). which is one of the most common light- 
curve fitting techniques, is adopted in this work. The output 
of SALT2 are B-band apparent brightness (mj), stretch fac- 
tor (xi) and colour (c) term. The distance modulus (hsn) 
is then calculated as 



Msjv = {tub — Mb) + axi — pc 



(1) 



where Ms is the standardised absolute SNe la magnitude in 
the B-Band for xi = c = 0, 01 describes the overall stretch law 
for the sample, and /3 is the colour law for the whole sample. 
These parameters are determined by minimising the scatter 
in the distance redshift relation or the Hubble residual (HR) 

HR = hsn{Mb,o., P) — fi C osmo{z, H , fi m , fU) (2) 
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Figure 1. Four examples of host galaxy spectra (black SEDs) at rest-frame wavelength and with the best-fit model templates from 
GANDALF/PPXF (red SEDs) overlaid. The top and lower panels arc galaxies hosting SNe la with large and small stretch factor values, 
respectively, as given by the labels. Left and right panels are high and low S/N spectra, with values given by the labels. The redshift for 
each object is also indicated by the labels. 



where Ho is the Hubble constant, SIm is the density param- 
eter for matter and Qa is the density parameter for the cos- 
mological constant. We aim at detecting systematic trends 
in the derived Hubble residuals and adopt Ho=65 km s _1 
Mpc -1 , S! m =0.3 and Q.a~0-7. The resulting light-curve pa- 
rameters are M=-30.11, a=0.12 and /3=2.86. 

For a negative Hubble residual, i.e. fi CO smo > MSiv, the 
cosmology places the SNe la at a greater distance than the 
peak luminosity predicts. This also means that the SNe la 
peak luminosity has not been fully corrected and it is too 
bright to match the cosmology. 



2.2 Host galaxy properties 

The aim of this work is to analyse the host galaxy prop- 
erties of SNe la from a spectroscopic point of view, par- 
ticularly using absorption line indices. For this purpose we 
need clean galaxy absorption spectra, free from contaminat- 
ing emission lines. Absorption line indices are measured on 
the clean spectra and analysed using single stellar popula- 
tion m odels . We use the me t hod d etailed in iThomas et al.l 
(|2010h and Ijohansson et alJ (|2012T l of which we provide a 
brief summary in the following section. 



2.2.1 Kinematics and Lick indices 

We utilise the fitting-code GANDALF (jSarzi et al.ll2006T l. which 
is based on on the penal i sed pi xel-fitting (PPXF) method 
of ICappellari fc Emselleml |2004) , to obtain clean absorp- 
tion spectra. A brief description of the code is given 
here, while the reader i s refer red to ISarzi et al l l|200rj ) and 
ICappellari fc Emselleml 12004 ) for details. GANDALF/PPXF si- 
multaneously fits stellar templates and emission line Gaus- 
sians to galaxy spectra. The result is a separation of emis- 
sion and absorption spectra. As stellar templates we adopt 
the single stellar population spec t ral en ergy distributions 
(SEDs) of lMaraston fc Strombackl i|201ll ). Besides the clean 
absorption spectra, other useful outputs of GANDALF/PPXF 



are stellar and gas kinematics, E(B-V) dust reddening and 
emission line fluxes/equivalent widths (EWs). 

Fig. [TJshows four examples of host galaxy spectra (black 
spectral energy distributions, SEDs) with varying stretch 
factor values (xi) and signal-to-noise ratios (S/N) per pixel 
(~1 A at 4500 A) in the g-band of the spectra, as given 
by the labels. The corresponding best-fit stellar templates 
from GANDALF/PPXF (red SEDs) are also shown. Prominent 
emission lines are visible in the galaxies hosting high stretch 
factor values. 

The Lick system consists of definitions of absorption 
line i ndices for 25 prominent absorption featu res in the op- 
tical (jWorthev et al.lll994l : iTraeer et al.lll998h . We measure 
the Lick indices on the clean absorption spectra, free from 
contaminating emi ssion lines, using the latest definitions of 
iTraeer et al.l (|l998l ). The resolution of the spectra are down- 
graded ^_to_tlw_lj£kjT12Sj^soluypn (that of the models, ~8- 
10 A. IWorthev fc Ottavianil Il997l) prior to measuring the 
indices. The index measurements are then corrected for ve- 
locity dispersion broadening, using the velocity dispersions 
found by GANDALF /PPXF together with the best fit stellar 
templates. The Lick indices are measured on the best fit stel- 
lar templates broadened to the Lick/IDS resolution, both be- 
fore and after further broadening with respect to the velocity 
dispersion measurements. The difference between these two 
measurements gives the correction factor which we apply to 
the Lick indices measured on the galaxy spectra. Using the 
error vectors provided with each SDSS spectrum we estimate 
Lick index errors through Monte Carlo simulations. 

It is important to note that these are all standard 
procedures. Hence, our derived Lick indices can be com- 
pared directly to other galaxy samples in the literature. 
This holds for samples with flux-calibrated spectra, as 
the Lick indices are sensitive to the instrumental con- 
figuration. When comparing Lick indices from different 
samples based on non-flux-calibrated spectra, the com- 
mon procedure is to use standard-stars observed with 
both instrumental configurations to deter mine possible off- 
sets/calibrations in the Lick indices (e.g. Ijohansson et alJ 
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bold ). In this work we fit the TMJ models to the ob- 
served Lick indices. The TMJ mode ls are based on the flux- 
calibrated stellar library M ILES (|Sanchez-Blazguez et alj 
120061 ; IJohansson et al.ll2010l ) . Hence the models are no longer 
tied to the non- flux-calibrated Lick/IDS system, making 
standard star-derived offsets unnecessary when using flux- 
calibrated galaxy spectra. Lick index offsets are crucial when 
adopting stellar population models based upon non-flux cal- 
ibrated stellar librar ies, such as the Lick/IDS library (e.g. 
IWorthev et al.|[l994l ). due to the sensitivity to the instru- 
mental configuration of the Lick indices as mentioned above. 

2.2.2 Stellar population parameters 

From the derived Lick indices we determine the stellar popu- 
lation parameters age and total metallicity and the element 
abundance ratios accessible through integrated light spec- 
troscopy of galaxies (O/Fe, Mg/Fe, C/Fe, N/Fe, Ca/Fe and 
Ti/Fe), using the iter ative method described in detail in 
IJohansson et al.l (|2012l ). This method is based on the TMJ 
stellar population models of absorption line indices. These 
are single stellar population models with variable element 
abundance ratios for the 25 Lick indices. Since several in- 
dices respond to variations of the same element abundances 
we have developed an iterative method. A x 2 -minimisation 
routine is used at each step to find the best-fit model. 

First we determine the traditional light-averaged stellar 
population parameters age, total metallicity, and a/Fe ratio 
from indices sensitive to these three parameters only. In the 
subsequent steps we add in turn particular sets of indices 
that are sensitive to the element the abundance of which 
we want to determine. In each step we re-run the x 2 -fitting 
code with a new set of models to derive the abundance of 
this element. This new set of models is a perturbation to the 
solution found for the base set and is constructed by keeping 
the stellar population parameters age, metallicity, and a/Fe 
fixed while modifying the abundance of the element under 
consideration. The derivations of element abundances are it- 
erated until the values do not further change. At the end of 
the sequence we re-determine the overall \ 2 using all indices 
together and re-derive the base parameters age, metallicity, 
and a/Fe for the new set of element ratios. The whole pro- 
cedure is iterated until the final % 2 stops improving by more 
than 1%. 

An enhanced a/Fe ratio is characterised by a depres- 
sion in Fe and reflects the ratio between total metallicity 
and iron abundance. Since O dominates the mass budget 
of total metallicity, the a/Fe parameter derived through 
our procedure can be reasonably interpreted as O/Fe. In 
IJohansson et alj (|2012T ) we therefore re-named the parame- 
ter a/Fe to O/Fe under the assumption that this ratio pro- 
vides an indirect measurement of oxygen abundance, i.e. 

[O/Fe] = [a/Fe] (3) 

2.2.3 Stellar mass 

We derive stellar masses (M*) from photometry accord- 
ing to standard SED-fitt ing as in iDaddi et al.l (|2005l ) 
and iMaraston" et all (|2006h via normalization of the SED. 
The SED-fitting is performed using the HyperZ code 
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Figure 2. The relationship between the error in stellar popula- 
tion parameters and quality of galaxy spectra, in terms of S/N 
per pixel (~1 A at 45 00 A) in the g-band (remake of Fig. 4 from 
IJohansson et alj|2012l) . Contours r epresent 3802 SDSS early-type 
galaxies from the MOSES sample llThomas et al.ll2010l) for error 
in age (upper panel), total metallicity (middle panel) and [O/Fc] 
(lower panel). Orange lines connect the mean errors in bins of 
g-band S/N. while green dashed lines show the linear behaviors 
below S/N=30. The distribution of the full sample of SNe la host 
galaxies (292 objects) is presented by the histogram, scaled by 
multiplying the binned histogram-values with 0.18/60. The black 
vertical line represents the cut at S/N=6, which defines our final 
sample. 

of iBolzonella. Miralles fe Pelldl (|2000l ) with SDSS spectro- 
scopic redshifts and extinction corrected ugriz magnitudes 
as input. Masses are derived using different s ets of stellar 
population models , i.e. fr om iMarastonl l|2005h (M05) and 
iBruzual fe Charlotl (|2003l ) (BC03). As is well known, the 
main difference between these models is the different treat- 
ment of the Thermally Pulsing Asymptotic Giant Branch 
(TP-AGB) contribution in the M05 models. The TP-AGB 
affects the luminosity at near-IR wavelengths for galax- 
ies with stellar population ages between ~ 0.2 and 2 Gyr 
|Marastonll2005l ). 

To obtain as robust total M* as possible, we ex- 
clude internal red dening in the fitting procedure, following 
IPforr et al.l (|2012i ). These authors find that the inclusion 
of galaxy internal reddening as a fit parameter generates 
masses that are often underestimated with respect to the 
true stellar mass. This result occurs because dust reddening 
combined with young model ages produce young dusty solu- 
tions as best-fit , with low ma ss-to-light ratios (the age-dust 
degeneracy, sec Rcnzini 2006). Furthermore, we exclude sin- 
gle burst templates which sometimes further increase the 
risk of getting underestimated stellar masses. 

To compare with literature values, which usually in- 
clude reddening, we further calculate stellar masses for the 
following options: M05 models with reddening (M05 re d), 
BC03 models without reddening and BC03 model with red- 
dening (BC03 red ). 

2.2.4 Gas-phase metallicity 

Gas-phase metallicities are derived from the EWs of 
the emission lines measured by GANDALF/PPXF (see Sec- 
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Figure 3. Distribution in redshift for the full sample of 292 
SNe la (open histogram) and the final sample of 84 SNe la 
(dash- filled histogram). The orange area represents the redshift 
range covered by the early-type galaxy sample used in Fig. [2] 
(represented by the contours) to show the relationship between 
stellar population parameter errors and quality of corresponding 
luminosity-weighted galaxy spectra. 



tion l2.2.Tj) . Several calibrations for the relationship between 
metallicity and emission-line ratios are available in the lit- 
erature. Relative metallicities generally agree between the 
various calibrations, while absolute metallicit i es fre quently 
disagree. The calibration of lKewlev fc Dopital |2002l . KD02) 
is commonly used in the literature, but it is based on the 
[7Vm6583]/[OT7A3727] ratio. The large difference in wave- 
length between the [NIIA6583] and [OIIA3727] lines make 
this ratio sensitive to dust reddening. Instead we use the 
03N2 index (log(([OIIIX5007]/H/3)/([NIIX658S\/Ha)) 
which is not sensitive to dust reddening due to the prox- 
imity of the emission lines in both the [OIIIA5007]/-ff / 3 and 
[NIIA6583]/Ha r atios. The Q3N2 i ndex was calibrated with 
12 + \og[0/H] in lPettini fc Pagell |2004l ) (PP04) 

(12 + log[O/H])pp 04 = 8.73 - 0.32 x 03iV2 (4) 

iKewlev fc Ellison! | 20081 ) re-calibrated the metallicity from 
PP04 on to the KD02 metallicity and found the following 
relationship 

(12 + \og[O/H]) Kn02 = 159.0567- 

54.18511 x (12 + log[O/#])pP04+ 
6.395364 x (12 + log[O/H]) FP0 4- 

0.2471693 x (12 + log [O/H]f PP0A (5) 

Using Eq. 1 4151 we determine metallicities for objects with 
detected emission lines. The presence of an active galactic 
nucleus (AGN) may strongly affect the derived metallici- 
ties since a stellar ionising radiation field is assumed in the 
comm only used metallicity calibrations (|Kewlev fc Ellison] 
120081 ). Hence we will therefore not consider gas metallici- 
ties for purely classified AGN hosts, while transition objects 
between star-formation and AGN are treated with caution. 



Table 1. Numbers of the host spectroscopy sample before and 
after the different cuts. The second column gives the number of 
spectroscopically confirmed SNe la and the third column both 
spectroscopically and photometrically identified SNe la. 



Selection 


Spec. Confirm 


Total" 


Full 


138 


292 


LC fitter cuts 


116 


196 


Separation cut 


105 


183 


Host spec, quality cuts 


63 


84 



"Total number used in the main analysis including both 
spectroscopically and photometrically identified SNe la 

2.3 Sample cuts 

The numbers of supernovae for the host spectroscopy sam- 
ple are given in Table [TJ where "Spec. Confirm" gives just 
the spectroscopically confirmed SNe la, for the different se- 
lection cuts described in this section. 

The individual objects of the host spectroscopy sample 
are presented in Table[2](the full version of this table is avail- 
able in Appendix A), where we provide SDSS-II Supernova 
Survey identification numbers, the SNe classification tech- 
nique used (Sp=Spectroscopically, Ph=Photometrically) 
and host galaxy coordinates and redshifts. We also pro- 
vide sample cut designation, where f=final host spec- 
troscopy sample selected through the cuts presented below 
(m=Master sample described in Section [3. 4. 2[) . 

2.3.1 SN la cuts 

To avoid contamination by low quality data, we cut the sam- 
ple to exclude the lowest quality SNe la observations. The 
accepted \ 2 °f the light-curve fit probability from SALT 
is constr ained and we remove everything with a reduced 
y 2 > 3 jLampeitl et alj l2010bh . Following iLampeitl et all 
|2010bf ). we retain SNe la with SALT2 parameters inside 
the ranges -4.5<xi<2.0 and -0.3<c<0.6, and remove any- 
thing with particularly large xi/c uncertainties. Outside 
these ranges ILampeitl et al.1 (|2010bl ) found the derived pa- 
rameters to be unreliable. These cuts of the fitted light-curve 
(LC) parameters reduce the sample to 196 objects ( "LC fit- 
ter cuts" in Table [TJ). 

It is possible that parameter gradients and local devia- 
tions from the integrated parameters will affect the results. 
Hence we discard SNe la that are separated by more than 
0.15' from the host galaxy, which is three times the SDSS 
fiber diameter (0.05'). This cut produces a sample of 183 
objects. 

2.3.2 Host spectroscopy cuts 

Fig.[2]shows the errors, represented by contours, in the stel- 
lar population parameters age (upper panel) , total metallic- 
ity (middle panel) and [O /Fe] (lower panel) as a function of 
the median S/N ratios per pixel (~1 A at 4500 A) in the g- 
band. This is the same format as Fig. 4 from ljohansson et alj 
(|2012T l; the errors are for 3802 SDSS early-type galaxies of 
the sampl e of Morphologically Selected Ellipticals in SDSS 
(MOSES, ISchawinski et al.ll2007l ; iThomas et alll2010r i, de- 
rived using the same method based on absorption line in- 
dices adopted in this work (see Section I2.2|l . Orange lines 
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Table 2. Specifications for the host spectroscopy sample. 



ID 


Class." 


Ra 


Dec 


z 


Sample' 1 






(deg) 


(deg) 






691 


Ph 


329.7300 


-0.4990 


0.1310 


m 


701 


Ph 


334.6050 


0.7976 


0.2060 


f 


717 


Ph 


353.6280 


0.7659 


0.1310 




722 


Sp 


0.7060 


0.7519 


0.0870 


f 


739 


Sp 


14.5960 


0.6794 


0.1080 


f 


762 


Sp 


15.5360 


-0.8797 


0.1920 


m 


774 


Sp 


25.4640 


-0.8767 


0.0940 


f 



"Supcrnovae classification, Sp=Spectroscopic, Ph=Photometric 
(sec Section [2{ 

^Sample selection, f=final host spectroscopy sample (see 
Section l2.3l l. m=Master sample (see Section 13,4.21 1 
(This table is available in its entirety in Appendix A. A portion 
is shown here for guidance regarding its form and content.) 

connect median errors in bins of g-band S/N, while green 
dashed lines show the linear behaviors below S/N=30. The 
errors on stellar population parameters increase almost lin- 
early with decreasing S/N over the S/N-range covered by 
the MOSES sub-sample. The quality, in terms of S/N in 
the g-band, of the 292 host galaxy spectra of our sample is 
illustrated by the histogram in the lower panel of Fig. [2] 

The distribution in redshift-space for the full SNe la 
sample is shown in Fig. [3] (open histogram) together with 
the redshift-range covered by the MOSES sample (orange 
area). Only a small fraction of the SNe la sample falls in 
and below the redshift-range of the MOSES sample. Thus 
the majority of the SNe la host galaxies have a lower S/N 
than the MOSES sample (as seen in Fig. [2}. It is reason- 
able to assume that the linear behaviour of the error-S/N 
relationships can be extrapolated to lower S/N. However, at 
the lowest values this may not hold such that the errors in- 
crease significantly. To avoid including data with too large 
uncertainties in the stellar population parameters, we cut 
the sample at a specific S/N value to balance the quantity 
and the quality of the sample. By changing the S/N-limit by 
one unit we find that the sample size is consistently reduced 
by ~15% up to S/N=6, while above this value the sample 
reduction increases to >25%. Hence, we only include host 
galaxies with S/N>6 in the g-band (vertical black line in 
the lower panel of Fig. [2}. We also discover one spectrum 
showing quasar features and two objects with telluric con- 
taminated Mgb absorption features. These quality cuts of 
the host galaxy spectra produce a final sample of 84 objects 
("Host spec, quality cuts" in Table [lj. The redshift distribu- 
tion of this sub-sample is shown by the dash-filled histogram 
in Fig. [3] The final host spectroscopy sample covers redshifts 
up to ~0.2. 



3 RESULTS 

Relationships between host galaxy parameters age, metal- 
licity, element ratios, velocity dispersion and stellar mass 
with SNe la properties stretch factor and Hubble residual 
are presented in this section. We do not find any signifi- 
cant trends between the host galaxy parameters and SALT2 
colour, hence such relations are not further discussed. 

For several of the studied relations we derive least- 
square fits. These are weighted using errors in the depend- 
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Figure 4. BPT-diagram for the 59 galaxies with detected emis- 
sion lines, i.e. AoN>1.5 for Ha, [OIII](A5007) and [Nil] (A6584). 
Blue points are star-forming galaxies and green points exhibit 
AGN activity only. The solid curved line is the theoretical star- 
formation limit from lKewIev et"aT | J200ll . KeOl) and the dashed 
curved line is the emp irical separation of AGN and purely star- 
forming galaxies from IKauffmann et al.l J2003bl . Ka03). Objects 
that fall between these lines are transition objects, hosting both 
star-formation and AGN activity. The sepa ration of LINER and 
Seyfert AGNs from lSchawinski et all j2007l . S07) is indicated by 
the solid straight line. 

able parameter. Estimated errors of the fit-parameters, de- 
rived using standard routines for least-square fitting, are 
used for computing the significance of the derived fits. When 
stated so, the least-square fits are sigma-clipped at a 2a 
level, i.e., we first fit a line using all data points, remove 
the outliers deviating more than 2a from the derived fit and 
finally redo the least-square fitting using the new sample. 

3.1 Emission line diagnostics 

Emission line diagnostics are measured on the emission- 
line spectra separated from the galaxy absorption spec- 
tra (see Section [2.2. ip . We use the emission line EWs and 
amplitude-to-noise (AoN) ratios to classify each galaxy as 
passive ly evolving, star-forming or hosting an AGN. Fol- 
lowing iBaldwin. Phillips fc Terlevichl (|l98ll ) (BPT) the ra- 
tios [OIII](A5007)/H,S and [NII](A6584)/Hq separate star- 
formation and AGN activity. The AoN for [OIII](A5007), 
[NII](A6584) and Ha sets the limit for emission line detec- 
tion. Due to the generally low S/N of the galaxy spectra 
(see Sectio nl2.3ll. we allow an AoN lim it of 1.5, compared to 
3.0 us ed in IKauffmann et al.l (|2003bl ) and ISchawinski et al.l 
|2007l ). Hence galaxies are classified as passively evolv- 
ing if they do not have an AoN>1.5 for [OIII](A5007), 
[NII](A6584) and Ha. With this limit we find 25 out of the 
84 galaxies in the final sample (see Section |2.3|) to be pas- 
sively evolving, i.e. 30 % of the SNe la in our sample occur 
in passively evolving galaxies. 

Fig. [4] shows the location of the emission-line detected 
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Figure 5. Comparison between stellar velocity dispersion and the different stellar masses derived in this work, i.e. M05 models without 
reddening (M05, panel a), M05 models with reddening (M05 ret j, panel b), BC03 m odels without reddening (BC03, panel c), BC03 
model with reddening (BC03 re( j, panel d) and the masses from lLampeitl et al.l ll2010bh (L10, panel e). The data points are colour-coded 
according to their emission line classifications from Section |3. II and given by the labels in the upper right panel, where passivc=passivcly 
evolving, AGN=AGN activity and SF=star-forming. Least-square fits to the cx-mass relationships are also included (solid lines) with 
the slopes given at the top of the panels. The standard deviations of the residuals to these fits are indicated by the labels. Star-forming 
galaxies with high L10 masses (lower left panel) are high-lighted as large blue diamonds in all panels, i.e. these symbols represent the 
same objects in all panels. 



galaxies in the BPT-diagram. The so lid, curved lin e is th e 
theoretical star-formation limit from iKewlev et ail (|200ll ). 
i.e. galaxies that fall below this line form stars. The dashed, 
curved line is the empirical separation of AGN and purely 
star- forming galaxies from iKauffmann et ail (|2003bl ). Ob- 
jects that fall between these lines are transition objects, 
hosting both star-formation and AGN activity. For this 
work we are interested in detecting star-form ation activity 
and c onsequently label everything below the IKewlev et al.l 
|200ll ) line as star- forming (blue points). The fraction of 
star-forming galaxies is 41 out of 84 (^50 %) and 18 galaxies 
are AGN (~20 %, green points). The sep a ration of LINER 
and Seyfert AGNs from ISchawinski et all (|2007l ) is also in- 
dicated by the solid straight line. Twelve of the 18 AGN 
labeled galaxies are LINERs and six are Seyfert AGNs. The 
colour coding (blue=SF, green=AGN) is kept throughout 



this paper, adding a red colour for passively evolving ob- 
jects. 



3.2 Velocity dispersion and stellar mass 

In this section, we assess the stellar mass derivations and 
highlight ingredients in the SED-fitting that are responsi- 
ble for systematic uncertainties in the derived masses. This 
exercise is important as we will later on use stellar velocity 
dispersion as a proxy for stellar mass. The evaluation of the 
stellar mass derivations is also important since stellar mass 
is extensively used in the literature when studying SNe la 
host galaxies. 

In Fig. [5] we compare the velocity dispersion measure- 
ments from GANDALF/PPXF (see Section l2~2~Tj) to the various 
photometric stellar masses (Section 12. 2. 3|) . e.g. M05 models 
without reddening (M05, panel a), M05 models with red- 
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Figure 6. Comparison between ages derived using absorption line indices (TMJ, see Section l2.2,21 and the derived ages from each of the 
four SED-fitting cases, e.g. M05 models without reddening (M05, upper left panel), M05 models with reddening (M05 re( j, upper right 
panel), BC03 models without reddening (BC03, lower left panel) and BC03 models with reddening (BC03 re( j, lower right panel). The 
one-to-one relation is included in each of the four panels (solid lines). Data points are colour labelled according to their emission line 
classifications from Section |3. II 



dening (M05 re( j, panel b), BC03 models without reddening 
(BC03, panel c) and BC03 model with reddening (BC03 re d, 
panel e) . A comparison betwee n the velocity dispersion mea- 
surements and the masses from lLampeitl et al, | <|2010bl . L10, 
panel e) is also included. For this sample we find 81 out of 84 
objects in common with our final sample. Data points are 
labeled according to the emission line classifications from 
Section 13.11 Objects from the L10 sample that are clearly 
offset from the general mass-sigma relation to higher masses 
are highlighted as large, open diamonds in each panel, i.e. 
these symbols represent the same objects in all panels. Least- 
square fits to the (j-mass relationships are also included in 
Fig. [S] (solid lines), with the slopes given at the top of the 
panels. The standard deviations of the residuals to these fits 
are indicated by the labels. 

IPforr et all (|2012h make a comprehensive study of the 
robustness of stellar population parameters that are derived 
from applying SED-fitting techniques to galaxy photometry. 
They use mock galaxies with known input stellar mass, age, 
etc. and study which SED-fit procedure is able to recover 
the input masses at best. They find that stellar masses are 
best reproduced when reddening is excluded in the SED- 
fitting procedure. In Fig. [5] we do indeed find the tightest 
correlation with velocity dispersion for the M05 case without 
reddening, judging from the standard deviations of the fit 
residuals. It can further be seen that the M05 case without 
reddening make the most pronounced distinction between 
star-forming galaxies (blue points) having lower masses than 
galaxies without emission lines, which is expected. 

By comparing panel a and panel c in Fig. [5] we can 
evaluate the effect of the different models used in the SED- 
fitting. The BC03 models produce higher masses compared 
to the M05 models, in particular for star forming galaxies 



(blue data points). This result is well-known and due to the 
stronger TP-AGB phase in the M05 models that contribute 
significantly to the luminosity of young stellar populations 
(see Section \2. 2. 3[) . The lower amount of near-IR light from 
the TP-AGB phase at intermediate-ages in the BC03 models 
is compensated by finding a best fit at older ages, when the 
RGB takes over in producing near-IR light. The older age 
of the best fit implies then a higher M*. In addition, the 
RGB takes over at older ages in the Padova tracks assumed 
in the BC03 models with respect to the M05 models, which 
contributes to the same effect (see Maraston et al. (2006) 
for details). 

The effect of including reddening in the fit can be ap- 
preciated by comparing panel a(c) to panel b(d) in Fig. 
In this case the inclusion of reddening results in lower ages, 
hen ce, lower masses in particular for star- forming galaxies 
(see lPforr et aljfeoij ). a significantly larger scatter in the a- 
mass relationship and consequently a flatter slope. The L10 
masses are derived using the Pegase models and including 
reddening in the fit. These masses show a relationship with 
velocity dispersion similar to the BC03 with reddening case, 
also because of the si milarity between the Pegase and BC03 
theoretical SEDs fsee lMarastonll2005t ). 

The objects with particularly high L10 masses com- 
pared to the velocity dispersion measurements, large open 
diamonds in each panel, shift significantly from panel to 
panel. The masses of these objects are significantly lower 
for the M05 without reddening case compared to L10. 

In Fig [6] the ages derived using absorption line indices 
(TMJ, see Section [2. 2. 2 1) are compared to those derived from 
each of the four SED-fitting cases, e.g. M05 models without 
reddening (M05, upper left panel), M05 models with red- 
dening (M05 re( j, upper right panel), BC03 models without 
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Figure 8. Relationship with host galaxy age (TMJ, panel a), velocity dispersion (panel b), stellar mass (M05 without reddening case, 
panel c), Z/H (panel d), O/Fe (panel e) and C/Fe (panel f) for the SALT2 stretch factor xi. The data points arc coloured according to 
the emission line classification of the host galaxies from Section 13.11 i.e. blue=SF, green=AGN and red=passivcly evolving. Solid lines 
are one time sigma-clippcd (2a level) least-square fits and the over-plotted crosses are removed data points. The slopes of the fits are 
given at the top of the panels. The value in parentheses in panel a is an additional error estimate for the slope (see text for detail). 
Horizontal lines are median xi-values in bins of the x-axis parameters with lengths indicating the width of the bins. Error bars in the 
lower left corners are average la errors. 
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Figure 7. Distribution of the stretch factor xi for the differ- 
ent emission classifications (sec Section 13.11 1 . The blue histogram 
represents star-forming, green AGN and red passively evolving 
galaxies. 



reddening (BC03, lower left panel) and BC03 models with 
reddening (BC03 re d, lower right panel). The one-to-one re- 
lation is shown in each of the four panels (solid lines). The 
best agreement with the TMJ ages are found for the M05 
case without reddening. This is also the case that shows the 



strongest correlation with velocity dispersion (see Fig. f5J). 
The effect of model choice (compare the upper left to the up- 
per right panel) is such that the BC03 models produce older 
ages for star-forming galaxies when compared to the TMJ 
ages. This result is due to the TP-AGB phase as already ex- 
plained. The inclusion of reddening in the SED-fitting lead 
to an underestimation of the ages (compare the left to the 
right panels) with respect to those derived from absorption 
line indices. 

Note that the agreement between TMJ ages and M05- 
SED-fit ages is remarkable and not simply due to the fact 
that the underlying population model is the same. Indeed, 
the results of SED-fit with reddening are in disagreement in 
spite of the use of the same model. This good result between 
photometricall y-derived and optical-absorption derived ages 
is not found bv lHansson. Lisker fc Grebel (|2012l ) in compar- 
ing ages from Lick indices from Gallazzi et al. (2005) based 
on BC03 models and SED-fit based on the same models. 



3.3 Stretch factor and stellar populations 

3.3.1 Final host spectroscopy sample 

Fig. [7] shows the distribution in stretch factor xi for star- 
forming (blue histogram), AGN (green histogram) and pas- 
sively evolving galaxies (red histogram), as classified ac- 
cording to their emission line strengths (see Section 13. 
It is clear that SNe la observed in passive galaxies have the 
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Figure 9. Left hand panel: binning in xi for stacking of spectra, where the limits of the bins are indicated by the horizontal dashed 
lines. Coloured labels state the number of objects within each individual bin. The histogram is distribution of xi. Right hand panel: 
resulting stacked spectra at rest-frame wavelength and for normalised flux density per unit A (f^), coloured according to the bin labels 
of the left hand panel. The labels, with corresponding colours, state the E(B-V) value for each stacked spectrum. 



lowest xi values hence the shortest decline rates, i.e. star- 
forming galaxies contain the most luminous SNe la. SNe la 
observed in galaxies with AGN activity appear to have xi 
values falling in between those with star-forming and pas- 
sively evolving host stellar populations. 

Relationships between the stellar populations parame- 
ters (TMJ age, Z/H, O/Fe and C/Fe), velocity dispersion 
and stellar mass with SALT2 stretch factor xi are presented 
in Fig. [8] The data points are colour coded according to 
the format in Fig. [7] The stellar masses are from the M05 
without reddening case, since they showed the best agree- 
ment with stellar velocity dispersion (see Section r3.2p . Solid 
lines are sigma-clipped least-square fits (see beginning of 
Section [3]) . The slope and corresponding errors of the fits 
are given at the top of the panels and sigma-clipped data 
points are indicated by over-plotted crosses. Horizontal lines 
are median values in bins of the x-axis parameters, where the 
length of the lines indicate the width of the bins. Mean la 
errors are shown by the error bars in the lower left corners. 

Clear anti-correlations between stretch factor and stel- 
lar population age (panel a, Fig. [H]), galaxy velocity disper- 
sion (panel b) and stellar mass (panel c) are found. Hence 
higher stretch factors, i.e. more luminous SNe la, are found 
in younger stellar populations and lower stellar masses. We 
find the slopes of the least-square fits to be different from 
zero at a >6a, 4<r and >3<r level for the xi-log(age), xi- 
mass and xi-log(cr) relationships, respectively. The fits are 
weighted with the xi-errors. However, in panel a the errors 
in age seem to dominate the scatter. To evaluate the esti- 
mated error on the fitted slope in the xi-log(age) relation we 
perturb the sample using errors in both xi and age through 
1000 Monte-Carlo simulations. For each realization we per- 
form a least-square fit and determine the standard deviation 
of the derived set of slopes. This value is given in parenthesis 
in panel a and is very close to the original error estimate. 



The stretch- factor also anti-correlates with Z/H (panel 
d, Fig. [5]) and C/Fe (panel f), but these relations are sig- 
nificantly weaker than those for age, mass and velocity dis- 
persion. The slopes of the least-square fits are shallower and 
differ from zero at less than a la level for Z/H and at a 1.5a 
level for C/Fe. The xi-O/Fe relation (panel e, Fig. [8]) is fiat 
with a large error on the fitted slope. 

To summarize, out of the parameters studied the stretch 
factor seems to primarily depend on the age of the host 
galaxy stellar populations. This age dependence is well in 
line with the correlation with star formation fraction dis- 
cussed above. 

Using the method described in Section \2. 2. 2 1 we also de- 
rive the stellar population parameters N/Fe, Mg/Fe, Ca/Fe 
and Ti/Fe. We find no strong correlations for xi with these 
parameters that would add valuable information to the anal- 
ysis. Measurements of N/Fe, Ca/Fe and Ti/Fe require higher 
S/N than the other parameters, while th e derived Mg/Fe 
and O /Fe ratios follow each other closely (|johansson et al] 
2012). Hence we do not further discuss these abundance ra- 
tios. 



3.3.2 Stacked spectra 

Due to the low S/N of the individual spectra used in the 
previous section, in the following analysis we perform con- 
sistency checks by using stacked spectra with relatively high 
S/N. We stack the spectra of the final sample in bins of xi- 
Hence, in practice we stack spectra of similar galaxy class 
because of the correlation between x\ and emission line clas- 
sification (see previous section). The spectra are first de- 
redshifted using the SDSS spectroscopic redshifts and lin- 
early interpolated to the same wavelength binning. We then 
normalise to the median flux density (fx) contained in the 
rest-frame wavelength range 5000-5500 A. The stacking is 
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Figure 10. Relationship with host galaxy age (panel a), velocity dispersion (panel b), Z/H (panel d), O/Fe (panel e) and C/Fe (panel f) 
for SALT2 xi for the stacked spectra. The data points arc coloured according to the colours of the binned spectra in Fig. [9] The la errors 
for the x-axis parameters are indicated by the horizontal error bars and bin widths are shown by the vertical bars. The xi value for each 
data point is the average xi value in the corresponding bin. Dashed lines are the least-square fits for the corresponding relationships of 
the individual spectra from Fig. [8] 



finally performed by ta king the median flux density value in 
each wavelength pixel |Lee et al.l l201oh. Using the median 
value is a safeguard against contaminated data, telluric con- 
tamination, features of individual spectra, and do es not bias 
agains t the very highest S/N spectra. Following iLee et al.l 
(|20ld ) we estimate the error in each pixel of the median 
stacked spectra with 

(3/N)s = y/Y^S/N)! (6) 

where (S/N) 3 is the S/N of the stacked spectrum and (S/N)i 
is the S/N of the individual spectra. 

The xi range covered by the final sample (see Sec- 
tion 12. 3p is divided into four bins, chosen to include over 
a dozen number of data points in each bin. The left hand 
panel of Fig. [9] shows the widths of the bins separated by 
the dashed horizontal lines, together with the correspond- 
ing number of objects in each bin. The histogram is the xi 
distribution for the final sample. The procedures applied to 
the individual spectra to compute the stellar population pa- 
rameters (see Section 12.2.11 and I2.2.2[) are then applied to 
the stacked spectra. 

The right hand panel of Fig. [5] shows the stacked spec- 
trum for each bin, following the colour coding of the left 
hand panel. The E(B-V) values produced by GANDALF/PPXF 
(see Section ^. 2. 1|) for each stacked spectrum are given by the 
labels with corresponding colours. The stacked spectra show 
a number of distinct features and emission lines that clearly 
change as a function of stretch factor. The largest stretch 
factors are found in host galaxies with the most pronounced 
star formation activity and the highest dust extinction. The 



host galaxies of faint supernovae with small stretch factors, 
instead, are emission line- free and dust-free, and show strong 
absorption features in their spectra. 

An analogue to Fig. [8] Fig. [TJTJ shows the relationships 
between stretch factor and stellar population properties and 
galaxy velocity dispersion for the stacked spectra. The data 
points in each panel are coloured according to Fig.|9]and the 
vertical bars indicate the size of the bins. The xi values are 
the mean within each bin. lcr error bars are shown for each 
data point in the horizontal direction. The lcr errors are 
similar for log(age), [Z/H], [O/Fe] and [C/Fe], i.e. smaller 
than 0.1 dex. The least-square fits to the final sample from 
Section r3.3.1l are included for comparison (dashed lines). The 
horizontal range covered in each panel has been truncated 
compared to Fig. [8] to better resemble the parameter range 
covered by the data of the stacked spectra. 

Lower xi values show older ages (panel a, Fig. I10[) . 
higher velocity dispersions (panel b) and higher total metal- 
licities (panel d), in agreement with the case of the individ- 
ual spectra (see Section T3. 3. l|l . However, only for age (panel 
a) we see this trend clearly for all data points adjacent in 
xi-space and significantly above the la error level. For total 
metallicity the trend with xi is diminished by the error bar 
overlap, which is due to the short range covered by this pa- 
rameter (~0.1 dex). A short parameter range is also found 
for [O/Fe] (~0.1 dex, panel e, Fig.[lO]) and [C/Fe] (~0.1 dex, 
panel f), resulting in a significant error bar overlap. We see 
no clear trends for these parameters. Hence the result of the 
individual spectra from Section 13.3.11 is reproduced for the 
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Figure 11. Relationship with stellar population age (panel a), velocity dispersion (panel b), [Z/H] (panel d), [O/Fc] (panel e) and 
[C/Fe] (panel f) for Hubble residual. The data points are coloured according to the emission line classification of the host galaxies from 
Scction l3.ll i.e. blue=SF, green=AGN and red=passively evolving. Solid lines are one time sigma-clipped (2a level) least-square fits and 
over-plotted crosses are removed data points. The slopes and corresponding errors are given at the top of the panels. Horizontal lines are 
median Hubble residuals in bins of the x-axis parameters with lengths according to the width of the bins. Error bars in the lower left 
corners are average 1<t errors. 



stacked spectra, i.e. xi show the strongest dependence on 
stellar population age. 



3.3.3 Comparison with the literature 

In agr ee ment with Oemler fc Tinslev (1 19791). van den Bergbl 
lll99d). iMannucci et all (|2005l ). ISullivan et all (|2006t ) and 
ISmith et alj (|2012T l, we find a higher fraction of SNe la 
events in star-forming compared to passively evolving galax- 
ies. Star-forming galaxies also show slower decline rates, 
i.e more luminous SNe la, comp ared to passive l y evo lving 
galaxies, a pattern fi rst noticed bylSullivan et al.l d2006t) an d 
further confirm e d bvlHowell et alj (120091) iNeill et all (12001) . 
iLampeitl etafl (|2010bl ) and lSmith et ail (|2012l ). For the first 
time we establish an anti-correlation between host galaxy ve- 
locity dispersion and stretch factor, suggesting faster decline 
rates in more massive galaxies. A similar anti-c orrelation is 
also f o und for stellar mass, in agre ement with iKellv et al.l 
(|20ld ). lLampeitl et al.l (|2010bl ) and lSullivan et al.l (|2QloT v 

For stellar metallicity and element abundance ratios we 
instead only find weak SALT2 xi dependencies. 

Furthermore, we find a clear anti-correlation between 
stretch factor and luminosity-weighted stellar population 
age which is more prominent and well defined than found 
for velocity dispersion and stellar mass. This result is true 
for both individual objects and stacked host galaxy spectra 
and indicates that the xi-mass relationship is a result of the 
corr elation between gala xy m ass and stellar populati on age. 

Hamuy et all (|2000l ) and I Gallagher et~al1 |200cJ ) inves- 



tigated absorption line indices allowing for distinction be- 
tween age and metallicity effects. The former prefer metal- 
licity over age as the main driver of S Ne la luminosity 
using only five objects. iGallagher et alj |20Qg| ) studied 29 
SNe la host galaxies and found age to be the dominant stel- 
lar population parameter affecting SNe la luminosi t y. Tw o 
recent studies, iHowell et ail |2009l ) and INeill et~aH (|2009l) . 
derived host galaxy parameters from photometry. Both au- 
thors favour age over metallicity as the SNe la luminosity- 
dependent factor, but determine the latter only indirectly by 
usin g derived masses and t he mass-metallicity relationship 
from Tremonti et al.l (|2004l ). In a recent study iGupta et ail 
|2011 ) also found faster decline rates for older stellar popu- 
lations using photometry, but do not include metallicity or 
element abundance ratios in their study. 

To conclude, this analysis strengthens the emerging 
trend in the literature that host stellar population age is 
the main driver of SNe la light-curve shape and luminosity, 
such that more luminous SNe la events occur in galaxies 
with younger stellar populations. 



3.4 Hubble residual and stellar populations 

3.4-1 Final host spectroscopy sample 

Fig.[TT]shows Hubble residual as a function of stellar popula- 
tion age (panel a), [Z/H] (panel d), [O/Fe] (panel e), [C/Fe] 
(panel f) and velocity dispersion (panel b) for the final host 
spectroscopy sample. The same colour coding and symbols 
as in Fig. [8] are used. We do not find significant trends for 
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Table 3. Stellar populations parameters for stacked spectra in bins of Hubble residual. The stacking limits are given in column 1, the 
number of objects in each stack in column 2 and the mean Hubble residual in column 3. The stellar population parameters together with 
corresponding errors are given in columns 4-7 and velocity dispersion in column 8. 



Stack 


N 
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log (age) 
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(Gyr) 
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43 


0.13 
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Figure 12. Relationship between gas-phase metallicity and Hub- 
ble residual for purely star-forming galaxies (blue solid points) 
and transition galaxies (open blue points). The solid line and 
the dashed line are least-square fits to all data points and to solid 
data points, respectively. The slopes and corresponding errors are 
given by the labels. Horizontal lines are median Hubble residuals 
in bins of the x-axis parameters with lengths according to the 
width of the bins for all data points. 



any of the parameters studied (<2a level for all parameters). 
To be sure that the quality of the host spectra does not di- 
lute possible trends, we have also stacked the spectra in two 
Hubble residual bins dividing the sample into sub-samples of 
roughly equal sizes. The stacking follows the procedure de- 
scribed in Section r3.3.2l Table|3]presents the Hubble residual 
range, number of objects and mean Hubble residual for the 
two bins, together with the stellar population parameters 
and velocity dispersion derived for the stacked spectra. We 
do not find significant differences for any of the stellar popu- 
lation parameters as well as for velocity dispersion between 
the two bins. 

The relationship between Hubble residual and gas- 
phase metallicity (see Section \2. 2. 4|) is presented in Fig 1121 
As pointed out in Section r2.2.4l AGN activity may affect the 
derived metallicities. The gas-phase metallicities are thus 
shown for galaxies with detected emission lines and classi- 
fied as star-forming. Galaxies in the transition region be- 
tween star formation and AGN activity (see Section T3.ip are 
plotted as open symbols. Least-square fits are derived for 
all data points (solid line) and purely star-forming galaxies 
(solid points, dashed line) with slopes given by the labels. 
Since this is a small sample, we do not sigma-clip the least- 
square fits. Median values in bins of metallicity are shown 
for all data points as horizontal lines. These lines show a 
clear trend of decreasing Hubble residuals for higher gas 
metallicities. The least-square fits also indicate a decrease 



in Hubble residual with gas-phase metallicity with a signif- 
icance of ~2.5a. 



This low statistical significance is caused by the two 
low-metallicity, low-Hubble residual data points. Consider- 
ing the individual Hubble residual errors both points deviate 
more than 2a from the least-square fits. Ignoring these two 
outliers in the least-square fitting results in a significantly 
steeper slope of -1.58±0.32. One possible explanation for the 
two outliers could be a failure in the emission-line template 
fitting (see Section 12.2.11) . affecting the derived gas- phase 
metallicities. However, this is not the case for either of the 
two objects. We have further studied the location of these 
SNe events with respect to the SDSS fiber location of the 
host galaxy. For the lowest-metallicity object the SNe event 
is within the SDSS aperture and it is a spectroscopically 
confirmed la. The other outlier object, not spectroscopically 
confirmed la, is located in the outskirts of the galaxy, at a 
distance of about two times the SDSS aperture diameter 
(3") from the fiber center. This could possibly explain the 
deviation of this object if a significant positive gas-phase 
metallicity gradient or a local metallicity upturn is present. 

Fig. 1 131 shows the relationship between Hubble residual 
and the stellar masses derived for the different prescriptions 
described in Section [2.2.31 e.g. M05 models without redden- 
ing (M05, panel a), M05 models with reddening (M05 re d, 
panel d), BC03 models without reddening (BC03, panel b) 
and BC03 model with reddening (BC03 re ,j, panel e). The 
relationship between Hubble residual and the L10 stellar 
masses is also included (panel f). The data points are colour 
coded, as usual, with blue=star-forming, green=AGN and 
red=passively evolving galaxies. Solid lines are least-square 
fits that have been sigma-clipped one time at a 2a level (see 
beginning of Section [3]). The slope of the fits are given at 
the top of the panels and sigma-clipped data points are in- 
dicated by over-plotted crosses. Horizontal lines are median 
values in bins of the x-axis parameters, where the length of 
the lines indicate the width of the bins. 

We can see in Fig. [T3] that the method of stellar mass 
derivation affects the Hubble residual-mass relationship. 
The masses derived with the M05 models without redden- 
ing show the least significant trend with Hubble residual 
and hence the smallest difference in Hubble residual for 
high and low mass objects. These masses also showed the 
strongest correlation with velocity dispersion in Section ^. 21 
In comparison, the masses derived with reddening and the 
BC03 models display stronger trends with Hubble residual 
but are weaker than for the L10 masses. The significance 
of the slopes of the least-square fits vary between 1.7a and 
2.2a. 
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Figure 13. Relationship between Hubble residual and the stellar masses derived for the different prescriptions described in Scction [2.2..3l 
e.g. M05 models without reddening (M05, panel a), M05 models with reddening ( M05 ret j, panel d) , BC03 models without reddening 
(BC03, panel b), BC03 model with reddening (BC03 re( i> panel e) and masses from iLampeitl et alj teOlObh (L10, panel f). The data 
points arc coloured according to the emission line classification of the host galaxies from Section 13.11 i.e. blue=SF, green=AGN and 
red=passivcly evolving. Solid lines are one time sigma-clippcd (2cr level) least-square fits and over-plotted crosses are removed data 
points. The slopes and corresponding errors are given at the top of the panels. Horizontal lines are median Hubble residuals in bins of 
the x-axis parameters with lengths according to the width of the bins. Error bars in the lower left corner arc average lcr errors. 



3.4-2 Extended photometric sample 

In addition to the SNe la with both available host galaxy 
SDSS spectroscopy and photometry, we can add SNe la with 
available SDSS photometry only to extend the sample for 
studying stellar masses. For this purpose we add the sam- 
ple studied in L10. These authors did a number of quality, 
redshift and host galaxy type cuts that produced a sam- 
ple of 162 objects. We followed the quality cuts of L10 for 
the host spectroscopy sample studied in this work, but we 
have also applied further quality cuts as described in Sec- 
tion [273] These additional constraints are lifted here in order 
to have a homogeneous sample. Specifically, the cut in S/N 
of the host galaxy spectra can be removed when we focus on 
photometry only. For the purpose of having a homogeneous 
sample we now also apply the redshift cut from L10, i.e., 
selecting objects with z<0.21. In fact, the S/N cut of the 
final host spectroscopy sample removed objects with z>0.21 
(see Fig. El). 

The host galaxy type cut in L10 was made in order to be 
able to study two distinct types of host galaxies, passive and 
star-forming, by constraining the allowed error on derived 
star-formation rates. Since we are not interested in distinct 
types of galaxies in this work, we can ignore this cut. The 
host spectroscopy sample used in our study is drawn from 
the same parent sample as the L10 sample, we will therefore 
mainly just add galaxies rejected due to the galaxy-type cut 
in L10. This action results in a sub-sample of 145 objects 
drawn from the full host spectroscopy sample. Out of this 



sub-sample 60 are in common with the L10 sample, leaving 
a Master photometric sample consisting of 247 objects. 

Fig. Q3] presents Hubble residual as a function of host 
stellar mass, derived with the M05 without reddening op- 
tion, for the Master sample. The sub-samples (with M05 
masses) are indicated by magenta points for the L10 sample 
and black points for the host spectroscopy sample with L10 
limits. Open, black-filled magenta points are overlapping ob- 
jects. The lines represent least-square fits to the different 
samples, i.e., black dashed line for the Master sample, black 
solid line for the host spectroscopy sample and magenta line 
for the L10 sample. The different samples were not individu- 
ally sigma-clipped in order not to bias the results by remov- 
ing different objects for the different samples. Instead, the 
Master sample was sigma-clipped twice at a 2a level and 
the objects removed through this step were ignored when 
performing the least-squares fitting of the two sub-samples. 
Green horizontal lines are median Hubble residual values for 
the Master sample, in bins of mass with lengths according 
to the width of the bins. 

It is clear that the slope of the fits, given by the labels, 
depends on the sample selection, where the L10 and host 
spectroscopy samples show the most and least prominent 
trends, respectively. The slope of the Master sample lies in 
between the values of the two sub-samples. The statistical 
significance of the slope is also sample dependent. The slope 
of the host spectroscopy sample is different from zero at a 
2.2<r level, while the corresponding numbers are 5.6a and 
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Figure 14. Relationship between Hubble residual and stellar mass for the different samples studied (sec text). The masses are derived 
using the M05 models without reddening (see Section 12.2,311 , The different samples are indicated by the upper labels with corresponding 
colours. , i.e. black points=host spectroscopy sample, magenta points=L10 sample and black filled magenta points=overlap sample. The 
lines show least-square fits to each sample, i.e. solid black line=host spectroscopy sample, magenta line=L10 sample and dashed black 
line=Master sample. The slopes of the fits and corresponding errors are given by the lower labels. Over-plotted crosses are data points 
removed by twice sigma-clipping the Master sample. Green horizontal lines are median Hubble residual values for the Master sample, in 
bins of mass with lengths according to the width of the bins. 



Table 4. Slopes and corresponding errors of least-square fits for the different samples studied and different stellar mass derivations. 



Sample 


M05 


M05 red 


BC03 


BC03 ret j 


Pegase (L10) 


Spec. 


-0.051±0.023 


-0.065±0.025 


-0.057±0.025 


-0.064±0.024 


-0.054±0.019 


L10 


-0.084±0.015 


-0.082±0.015 


-0.072±0.013 


-0.076±0.013 


-0.077±0.014 


Master 


-0.064±0.012 


-0.067±0.013 


-0.062±0.013 


-0.071±0.012 


-0.068±0.011 



5.3cr for the L10 and Master sample, respe ctively. For the 
L10 sa mple we recover the trend found in iLampeitl et al.l 
(|2010bD , but with a slightly steeper slope due to the sigma- 
clipping procedure. We have repeated this analysis and con- 
firm the above results for the different stellar mass deriva- 
tions, including the procedure used in L10. Table [4] presents 
slopes and errors of least-square fits for the different samples 
and different stellar mass derivation. The systematic varia- 
tions apparent for the smaller sample (see Section [3ATJ are 
weaker when the number of galaxies is increased. It should 
also be noted that the Hubble residuals for the L10 sample 
have been re-derived using the a and ft values used in this 
work (see Section (2T|. 

Two effects are responsible for the variation in slope of 
the different samples; 1. The median Hubble residual val- 
ues (green horizontal lines in Fig. 1 14 p show a drop around 



log M0 = 1O.2 and are close to constant below and above this 
value. 2. The L10 and host spectroscopy samples have differ- 
ent stellar mass distributions, which is apparent in Fig. 1141 
This is better emphasized in Fig. 1151 where the top panel 
shows the stellar mass distributions for the different sam- 
ples, i.e., magenta and black histograms represent the L10 
and host spectroscopy samples, respectively. The host spec- 
troscopy sample is lacking low-mass objects, while the L10 
sample is deficient in high-mass galaxies. The high-mass bias 
of the host spectroscopy sample is not surprising as SDSS 
spectroscopic targeting r equired a magnitud e of r<17.7 for 
the main galaxy sample l)Strauss et al.ll2002h . 

The discrepancy in slope between the different samples 
is driven by the different mass ranges covered. This result 
is shown in the main panel of Fig. 1151 which is a remake 
of Fig. 1141 Here we instead study the Hubble residual dis- 



© 2011 RAS, MNRAS 000, [[>?? 



SNe la SDSS-II host galaxy properties 17 




9 10 11 

M05 Stellar Mass [log(M )] 



Figure 15. Upper panel: Histograms of the stellar mass distribution for the different samples studied (see text) as indicated by the 
labels. Main panel: Re-make of Fig, 1141 but also adding high redshift objects (z>0.21, cyan points, not used in the least-square fits). 
Vertical dashed lines indicate boundaries between different mass-ranges for which we derive least-square fits. The slope of each fit is 
given by the labels with corresponding colours. Over-plotted crosses are removed data points through one-time sigma-clipping. Right 
panel: The Hubble residual distribution for each stellar mass bin of the main panel. The histograms are colour-coded according to the 
top labels. Lower labels give median values in each mass bin for the same colour coding. 



tributions in three different mass intervals: log M© <9.5, 
9.5< log M© <10.2 and log M© >10.2. Least-square fits, de- 
rived without the cyan points, in these different regions are 
presented by the lines with slopes according to the labels 
with corresponding colours. For the high-mass bin the HR- 
mass relation is clearly flat. This appears to be the case also 
for the low-mass bin, but this regime is more sparsely pop- 
ulated. Higher Hubble residuals are found at low masses, 
while in the region around logM© ~10.0 the Hubble resid- 
ual varies significantly. Hence, the Hubble residual-mass re- 
lationship shows a step function where negative and positive 
values are found for high- and low-mass objects, respectively. 
This behaviou r can also be found in the sample studied in 
ISullivan et all (|2010l ). R is clear from their Fig. 4 that the 
Hubble residuals suddenly drop above logM© ~10.0, while 
above and below this limit the average Hubble residuals stay 
roughly constant. Stellar masses around 10 9 M© may not 
represent the total mass of a galaxy, but the mass of the 
star-forming fraction, which outshines the underlying older 
popu lation (the "tip of th e iceberg" effect, Maraston et all 
|2010| ). The simulations of IPforr et all (|2012h indicate that 



when a lO n M0 galaxy has experienced a 1-10 % by mass 
recent starburst, the SED-fit derived mass is closer to the 
mass of the burst rather than to the total mass. This is an 
interesting effect to keep in mind when a large mass range 
is considered. 

We test the significance of the step function as com- 
pared to the linear fit over the full mass-range using the 
resulting x 2 -values and an F-test for nested regression mod- 
els. To have the same sample for both models we use the 
full Master sample (247 objects) and take into account that 
the step function has six parameters instead of two for the 
linear fit. We find that the F-test rejects the null hypothesis 
that the step function does not provide a significantly better 
fit than the linear fit at a probability >90% (93.5%). 

The Hubble residual distributions of the different mass 
bins are plotted in the right panel of Fig. \TE\ with colours 
according to the labels at the top of this panel. The me- 
dian values of each mass bin is given by the lower right 
labels. The median Hubble residual of the low-mass bin is 
offset from the high-mass bin by 0.147 mag, while the me- 
dian value of the intermediate mass bin is located between 
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these values. A Kolmogorov-Smirnov test confirms that the 
distributions of the high- and low-mass bins are different at 
a significance level of over 99.9%. Hence, the stellar mass dis- 
tribution of the selected sample strongly affects the derived 
Hubble residual-mass trend. This would not be the case if 
there was a continuous trend over the full stellar mass range. 

To further emphasize this point we added the sub- 
sample of objects with z>0.21 (cyan points) in Fig. 1151 In the 
top panel it is clear that this sample is distributed towards 
high stellar masses. This sub-sample exhibits low Hubble 
residuals similar to the high-mass bin of the Master sample, 
as can be seen in the main and right panels. Including this 
sub-sample in the Master sample significantly reduces the 
slope of the least-square fit from -0.065 to -0.048. 

We confirm similar step functions for all stellar mass 
derivations, i.e. M05 with reddening, BC03 without redden- 
ing, BC03 with reddening and also for the stellar masses 
derived with the setup of L10. We have also applied the 
cuts described in Section \2. 31 to the Master sample and L10 
sample and find that the results are not affected. However, 
since the cut in spectral resolution can not be applied to the 
L10 sample, it has not been considered for this consistency 
check. 



3.4-3 Comparison with the literature 

We find that Hubble residual depends on host galaxy mass 
in a greement with several recent reports in the litera- 
ture (IKellv et al.ll2010l;lLampeitl et al.ll2010bl ; ISullivan et all 
l2010l ; lGupta et alj|201ll ). However, we find that the Hubble 
residual-mass relationship beh aves as a step function , which 
can also be seen in the work of ISullivan et al.1 (^Old ). where 
the Hubble residuals mainly change around log M ~10.0. 
For the first time, we study the relation between Hubble 
residual and stellar velocity dispersion, but find no signifi- 
cant trend. However, this may be due to the high-mass bias 
of the sample studied. 

It is believed that the Hubble residual-mass trend is the 
result of a more fundamental relationship between Hubble 
residual and SNe la progenitor properties that could pos- 
sibly be mimicked through host stellar populations. Sev- 
eral authors have therefore also studied relationships be- 
tween Hubble residual an d the stellar po p ulatio n param- 
eters age and metallicity. iGallagher et al.l (|2008l ) presents 
the only study in the literature to date that is based on 
absorption line indices and find stellar metallicity to be 
the source of systematic Hubble residual variations for a 
small sample o f ~20 host galaxies. I Howell et al.l (|2009T ) and 
ISullivan et ail (|2010T l instead infer gas-phase metallicities 
from stellar ma sses together w it h the mass-metallicity re- 
lationship from iTremonti et al.l |2004h and find significant 
trends in these parameters with Hubble residual. Gas-phase 
metalli cities are also infe rred from stellar mass in a similar 
way bv lNeill et al.l (|2009j ), but they instead find a stronger 
relation between ste llar population age and Hubble residual. 
Similar to this work. lD' Andrea et al.l ( 201 ll ) derive gas-phase 
metallicities directly from the emission line ratios and find a 
trend of higher Hubble residuals in galaxies of lower metal- 
licity. Hayden et al. (in prep.) find that the scatter in the 
Hubble residuals are significantly reduced when gas-phase 



metallici ty is added as furt her fit parameter besides stel- 
lar mass. iGupta et al.l ()201 lh instead find stellar population 



age to correlate with Hubble residual. For the sample of 84 
host galaxies studied in this work, we do not find significant 
trends of Hubble residuals with any of the absorption line 
derived stellar population parameters age, metallicity and 
element ratios. Again, this may be due to the high-mass 
bias of the sample studied. 

We have shown in the previous two sub-sections that 
the Hubble residual-mass trend is dependent on the mass- 
range spanned by the galaxies. The lack of a trend between 
Hubble residual and velocity dispersion as well as with stel- 
lar population age and metallicity could therefore be the re- 
sult of the limited mass-range studied through the final host 
spectroscopy sample. The reason why we find a tentative re- 
lationship for gas-phase metallicity probably arises from the 
fact the star-forming galaxies, for which we can infer gas- 
phase metallicities, cover the stellar mass range where the 
Hubble residuals dramatically vary (around 10 10 M©; com- 
pare top left panel of Fig. [T3] and Fig |14[l together with the 
fact that emission lines are more robustly derived in low-S/N 
spectra than absorption lines. 



4 DISCUSSION 

We derive host galaxy stellar population parameters for 
SNe la from the SDSS-II Supernova survey. These parame- 
ters are derived from absorption line indices through com- 
parison with stellar population model predictions (see Sec- 
tion Due to the low S/N of the host galaxy spectra we 
focus on the highest quality objects, considering both the 
SNe la observations and host galaxies, which results in a 
final sample of 84 objects (see Section f2 . 3 [) . 

We find a strong relationship between SALT2 stretch 
factor and stellar population age (see Section 13. 3p . For 
SALT2 colour on the other hand, we do not find dependen- 
cies on any of the host galaxy parameters, in agreement with 
previous studies. We also derive photometric stellar masses 
for the host galaxies using a variety of methods and models. 
The Hubble residual-stellar mass relation is found to behave 
as a step function, where the trend is flat at the high-mass 
end (see Section ^. 4. 2|l , This is also the mass-regime covered 
by our host spectroscopy sample, which may explain a lack 
of trends between Hubble residual and the stellar population 
parameters studied. 



4.1 Progenitor systems 

Considering several different progenitor systems, theo- 
retical models of SNe la explosions find SNe la rates 
SNRs) peaking at delay times below or close to 1 Gy r 
Yungelson fc Livioll2000l : iGreggiol [20051 ; iRuiter et al.ll2010l ) . 
The SNRs then smoothly decline and become 10-100 times 
lower at delay times of ~10 Gyr. Thus theory implies that 
the delay times of SNe la span the range from relatively 
quick events to explosions delayed by a Hubble time. 

Using the host galaxy stellar population ages as delay- 
time proxies, observations support the theory as we find a 
wide range of lumin osity-weighted ages from below 1 Gyr to 
>10 Gyr. Similarly [Gupta et al.l (120111 ) find SNe la events 
in stellar populations w ith mass- weig h ted a ges > 10 Gyr and 
down to ~2 Gyr, while iHowell et all (|2009f ) and lNeill et~al] 
(|2009f ) present luminosity- weighted ages from ~100 Myr up 
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to >10 Gyr. The work of lBrandt et all (|201Ch suggests delay- 
times of <400 Myr and >2.4 Gyr, indicative of two differ- 
ent progenitor channels. All these studies agree with the- 
ory such that a higher fraction of SNe la events are found 
in young star-forming galaxies compared to old passively 
evolving galaxies. If SN la peak luminosity increases with 
progenitor mass, the wider range of stellar masses in star- 
forming galaxies together with short life-times for massive 
stars result in the higher SNe la fractions. 

Indeed, considering the strong anti-correlation between 
stretch factor and stellar population age (see Section 13.31) , 
there is a clear connection such that those SNe la with the 
shortest delay times, hence most massive pr ogenitors, are 
the m ost luminous and vice versa. Contrary to lBrandt et al.l 
|20ld ). the relation between stretch factor and stellar pop- 
ulation age we found shows a smooth behaviour which may 
be i ndicative of a single proge nitor system. 

lYungelson fc Livid (|200d ) show that the SNR for a DD 
system peaks at delay times of ~100 Myr, while the analo- 
gous for a SD system is ~1 Gyr. Thus the lower age limit 
of the stellar populations hosting SNe la could be used to 
constrain possible progenitor systems. However, it may be 
unreasonable to attribute the integrated light of a galaxy 
to be a good proxy for SNe la delay times, considering the 
range of stellar populations, age and metallicity, possibly 
present within an instrumental aperture. An alternative ap- 
proach could be to use observations confined to the vicinity 
of the SNe la. However, it is likely that a mixing of stellar 
populations occur over a long time span. The probability of 
identifying the "true" parent stellar population for a SNe la 
progenitor should increase for younger stellar populations. 
It may therefore be possible to constrain the lower delay 
time limit using the stellar populations of the vicinity of the 
SNe la. This information can then in turn constrain possi- 
ble progenitor systems when comparing to theoretical SNe la 
models. 

4.2 Hubble residual minimisation 

The re ports of Hubble residual dependencies on host galaxy 
mass (Kelly et al.l l2010l: ISullivan et ail 120101 : iGupta et all 
|2011| ; IConlev et alj|201ll ) have led to the inclusion of this 
parameter to further reduce the scatter in the redshift- 
distance relation (e.g. iLampeitl et al.|[2010bh . In this work 
we find systematic uncertainties in the derived masses, due 
to the adopted SED-fitting prescription, to affect the Hub- 
ble residual-mass relationship (see Section [3.4. 1 [I . However, 
these systematic uncertainties become weaker when the 
number of galaxies is increased (see Section VS. 4. 2 [) , 

We also find that the stellar mass range covered affects 
the derived Hubble residual-mass trend (see Section [3.4.2p . 
This trend behaves as a step function, where the low- and, 
in particular, the high-mass end show flat trends. Due to 
this step function, the stellar mass range sampled will result 
in a systematic variation of the derived Hubble residual- 
mass slope. Hence, the non-linear behaviour of the Hubble 
residual-mass relationship should be accounted for when us- 
ing stellar mass as a further parameter for minimising the 
Hubble residuals. 

The dependence of Hubble residual on stellar mass is 
believed to be a consequence of a dependency on a more fun- 
damental parameter. The sampling of the high-mass regime 



of our host spectroscopy sample, where the trend between 
Hubble residual and stellar mass is flat, hampers a conclu- 
sion on which parameter is driving Hubble residual varia- 
tions. 

Interestingly, the step in the Hubble residual-mass plane 
appears at around ~2x10 10 Mq which is close to the evo- 
lutionary _transitk2ri_niass of low- redshift galaxies, discov- 
ered by iKauffmann et all (|2003al ). They find that below 
~3x10 10 Mq, galaxies are typically disc- like objects with 
young stellar populations and low surface mass densities. 
Above this limit, they instead find an increasing fraction 
of bulge-dominated galaxies with old stellar populations. 
Hence, the transition mass marks a change in the stellar 
populations of galaxies. The similar mass-range shared be- 
tween the sudden shift in the Hubble residual and transition 
in galaxy properties is therefore most likely not coinciden- 
tal. The stellar populations below and above the transition 
mass bare witness of the properties of the SNe la progeni- 
tors with high and low Hubble residuals, respectively. If the 
Hubble residual-mass relation is indeed flat at the low-mass 
end (only tentative due to low number statistics) as it is 
at the high-mass end, the step function may indicate two 
samples of SNe la with high and low Hubble residuals, i.e. 
showing a weak difference in peak luminosity not accounted 
for by the light-curve fitting techniques. 

The overlap in mass for the galaxy evolutionary transi- 
tion and change in Hubble residual tells us that the latter 
does not change significantly in old, passively evolving stel- 
lar populations. Instead there is a strong variation in lower 
mass galaxies with younger stellar populations, higher S FRs 
and lower gas-phase metallicities (|Tremonti et a,l.ll2004h . In- 
deed several authors suggest gas-phase metallicity to be the 
fundamental param e ter driving Hubble r e sidual variations 
(iHowell et all 120091 : ISullivan et all l20ld ; | D'Andrea et al.l 
l201ll . Hayden et al. in prep.). Also in this work we find tenta- 
tive results for such a relation between Hubble residual and 
gas-phase metallicity (see Section 13.4. 1|) . If this is indeed 
the fundamental parameter, it represents progenitor metal- 
licity. This is in agreement with theore tical predictions of 
prompt SNe la in star-forming galaxies |Yungelson fc Livid 
l200d ; iGreggiol 120051 ; iRuiter et all l20ld ). The flat slope at 
masses above ~2xlO lo M0 is then expected as the gas- 
hase metallicities sho w a flat behaviour in this mass-regime 
Tremonti et~aH 12004 Hayden et al. in prep.). However, 
prompt SNe la are not expected in old stellar populations 
present in passively evolving, massive galaxies. The deriva- 
tion of the full range of stellar population parameters and 
gas-phase properties for a high quality data set, covering 
the full mass range, is required to determine the fundamen- 
tal parameter driving Hubble residual variations. 

Future data from IFU spectroscopy, sampling the local 
environment of SNe la will be quite useful in identifying the 
fundamental parameter responsible for the scatter in the 
Hubble residuals. For such data sets it will be important to 
sample the full stellar mass range, in particular the mass 
range around 1O 1O M . 



5 CONCLUSIONS 

We present an analysis of the stellar populations of SNe la 
host galaxies using SDSS-II spectroscopy. Using the stel- 
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lar population mod els of absorption line indices from 
iThomas et al.l (|201 lh we derive the stellar population pa- 
rameters age, metallicity and element abundance ratios to 
study relationships with SNe la properties. We also mea- 
sure stellar velocity dispersion from stellar template fitting 
and gas-phase metallicity from emission lines when detected. 
Furthermore, we derive stellar masses and revisit the corre- 
lation between Hubble residual and stellar mass that has 
recently received much attention in the literature. The stel- 
lar masses are derived from SDSS-II photometry from repeat 
Stripe 82 observations and several SED-fit methods. From a 
sample of 292 SNe la we select 84 objects depending on the 
quality of the host galaxy spectroscopy and accuracy of the 
SNe la properties (see Section [2. 3 II . 

We find a larger fraction of SNe la together with typi- 
cally higher SALT2 stretch factor values (i.e. more luminous 
SNe la with slower declining light-curves) in star-forming 
compared to passively evolving galaxies (see Section T3.3p . in 
agreement with the literature. Hence, previous studies sug- 
gest that the decline rate and peak-luminosity of SNe la 
depend on stellar population age. With the large parame- 
ter space covered in this work, we indeed find that SALT2 
stretch factor values show a strong dependence on stellar 
population age in the sense of a clear anti-correlation (see 
Section I3.3[) . As a result of this trend we also find anti- 
correlations with velocity dispersion and galaxy mass, while 
only weak anti-correlations are also found for stellar metal- 
licity and the element abundance ratios studied. 

To ensure that the quality of the host galaxies are not 
affecting the selected sample, we stacked the spectra in bins 
of stretch factor. This exercise confirms the above results 
that stellar age is a strong candidate for the main driver of 
SNe la luminosity (see Section [3.3.2p . Hence, SNe la peak- 
luminosity is closely related to the age of the stellar progen- 
itor systems, where more luminous SNe la appear in young 
stellar populations. 

The peak-luminosity variation of SNe la is corrected for 
through light-curve fitting. However, scatter in the redshift- 
distance relation even after these corrections may intro- 
duce uncertainties in derived cosmological parameters. It 
has been debated which parameter that could be used to 
further minimize this scatter. While stellar mass has already 
been introduced for this purpose, stellar population age and 
metallicity as well as gas-phase metallicity have been pro- 
posed as fundamental parameters for the Hubble residual 
dependency. 

However, we identify no statistically significant trends 
of Hubble residuals with any of the stellar population pa- 
rameters studied or with stellar velocity dispersion (see Sec- 
tion 13.4. 1 f) . Instead we find tentative results confirming a 
trend between Hubble residual and gas-phase metallicity as 
previously reported. For the Hubble residual-stellar mass 
relationship, our selected sample shows a weak trend that 
is affected by the method of stellar mass derivation. Stel- 
lar masses that show the tightest correlation with stellar 
velocity dispersion produce the weakest trend with Hubble 
residual. To study this in more detail we extend the sam- 
ple to include 102 SNe la with available SDSS host galaxy 
photometry and lift the spectroscopic quality constraints, re- 
sulting in a sample of 247 objects for deriving stellar masses 
from SDSS photometry. With the better statistics of this 



larger sample the systematic variations arising from the stel- 
lar mass derivation are weakened (see Section ^. 4, 2[) . 

For the extended photometric sample it is clear that 
the reported Hubble residual-mass relation is strongly de- 
pendent on the stellar mass range studied and behaves as 
a step function. In the high-mass regime, the relation be- 
tween Hubble residual and stellar mass is flat. Since our sam- 
ple with available host galaxy spectroscopy mainly probes 
this high-mass regime, it is not surprising that we do not 
find any significant Hubble residual trend with the stel- 
lar population parameters studied. Below a stellar mass of 
~2x 1O 1O M0, i.e., close to the evolutionary transition mass 
of low-redshift galaxies reported in the literature, the trend 
changes dramatically such that lower-mass galaxies exhibit 
fainter SNe la after light-curve corrections. We conclude that 
the non-linear behaviour of the Hubble residual-mass rela- 
tionship should be accounted for if stellar mass is to be used 
as a further parameter for minimising the Hubble residuals. 
However, it is crucial to find the fundamental parameter 
driving systematic variations in the Hubble residual. 
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APPENDIX A: 

Table Al. Properties of the host galaxies 
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-0.8348 


0.1510 


f 


1748 


Ph 


353 1 1 20 


-0.4822 


0.3400 




1825 


Ph 


316.8200 


0.3590 


0.1560 


in 


1906 


Ph 


309 6490 

- ty ' •.' . V7 ^ (7 V 7 


0.9133 


2000 


[' 


2056 


Ph 


320.0370 


-0.3899 


0.1900 




2092 


Ph 


336 4560 

* 7 ' 7 ' V . i " 7 V7 V7 


-0.2390 


0.1430 




2331 


Ph 


16.6430 


1.0873 


0.2630 




2361 


Ph 


7.6140 


-0.4533 


0.1540 




2561 


Sd 

1 J F 


46.3440 


0.8597 


0.1190 


m 


2689 


Si) 
o f 


24 9000 

ijrt. CV7V7V7 


-0.7579 


0.1620 


f 


2766 


Ph 


52.7120 


-0.6186 


0.1500 


f 


2992 


Sd 
°F 


55 4970 

"7 : 7 ■ i 17 J V7 


-0.7829 


0.1270 




3241 


Sd 

1 'F 


312.6530 


-0.3542 


0.2590 




3293 


Ph 


39 6390 

' ) v ' . V7 (7 -7 \7 


0.2862 


0.1320 




3565 


Ph 


2.5900 


0.7100 


0.2900 




3592 


Sd 
°F 


19 0530 


0.7905 


0.0870 


[' 


3888 


Ph 


7.4840 


-0.4557 


0.1180 


[' 


3892 


Ph 


18.6860 


-0.4518 


0.3500 




3901 


Sd 


14.8500 


0.0026 


0.0630 


f 


4019 


Ph 


1.2620 


1.1464 


0.1810 


in 


4065 


Ph 


47.9980 


1.0513 


0.1310 


f 


4181 


Ph 


37.8160 


-1.1343 


0.3420 




4651 


Ph 


37.3760 


-0.7475 


0.1520 




4690 


Ph 


32.9300 


0.6877 


0.1990 


in 


4966 


Ph 


31 3 44QD 


-1.0414 


0.3130 




5230 


Ph 


334.7930 


0.8388 


0.3220 




5736 


Sp 


22.8660 


-0.6288 


0.1440 


in 


5785 


Ph 


328.5980 


0.0844 


0.1480 


f 


5944 


Sp 


29.2020 


-0.2126 


0.0460 


in 


5966 


Sp 


16.1910 


0.5133 


0.3100 




6057 


Sp 


52.5540 


-0.9745 


0.0670 


f 


6213 


Ph 


344.1060 


-0.4500 


0.1090 




6295 


Sp 


23.6740 


-0.6042 


0.0800 




6332 


Ph 


325.9660 


-0.7183 


0.1520 




6406 


Sp 


46.0890 


-1.0631 


0.1250 


in 


6558 


Sp 


21.7020 


-1.2382 


0.0570 


f 


6638 


Ph 


45.0390 


-1.2376 


0.3260 




6683 


Ph 


327.2140 


0.6126 


0.5090 





"Supernovae classification, Sp=Spectroscopic, Ph=Photomctric 
(see Section [2]l 

''Sample selection, f=final host spectroscopy sample (see 
Section l2.3l l. m=Master sample (see Section |3. 4. 21 1 
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Table Al - continued 



SN ID 


Class. a 


R,a 


Dec 


H.easniit 


Sampl 






( dc g) 


(dcg) 






6813 


Ph 


27.2700 


0.0566 


0.4110 




6851 


Ph 


52.1050 


-0.0488 


0.3050 




6962 


Sp 


38.8610 


1.0745 


0.0940 


f 


7335 


Sp 


318.8830 


-0.3546 


0.1970 


in 


7350 


Ph 


7.5630 


-0.7872 


0.1550 




7431 


Ph 


340.9550 


-0.2746 


0.3500 




7876 


Sp 


19.1830 


0.7936 


0.0760 


f 


8004 


Pli 


347.5290 


-0.5580 


0.3510 




8195 


Ph 


o o "i f\f\r'f\ 

331.0060 


-0.8957 


0.2690 




8280 


Ph 


8.5730 


0.7955 


0.1850 


in 


8555 


Ph 


2.9160 


-0.4150 


0.1980 


in 


8888 


Pli 


5.1660 


0.3214 


0.3990 




9117 


Ph 


A t ' f\f\f\f\ 

46.9000 


0.9884 


0.2720 




9133 


Ph 


16.6430 


0.4606 


0.2670 




9558 


Ph 


15.9840 


0.3810 


0.3910 




9633 


Ph 


33.9190 


1.0952 


0.1960 


in 


9739 


Pli 


323.6950 


-0.8790 


0.1200 




9817 


Pli 


5.0380 


0.5950 


0.2250 




10028 


Sp 


17.7420 


0.2762 


0.0650 


f 


10096 


Sp 


29.4300 


-0.1794 


0.0780 


f 


10434 


Sp 


329.9570 


-1.1924 


0.1040 


f 


10690 


I'll 


347.5000 


1.0822 


0.3120 




10805 


Sp 


344.9280 


-0.0136 


0.0450 


f 


11102 


I'll 


315.1590 


-0.9996 


0.1540 


f 


11120 


I'll 


325.4590 


0.3693 


0.1070 


f 


11172 


I'll 


322.4130 


-0.2022 


0.1360 


f 


11294 


I'll 


20.2970 


0.3737 


0.1170 


in 


11296 


Ph 


20.5190 


0.3106 


0.1760 


in 


11306 


Ph 


56.7380 


-0.5176 


0.2740 




11311 


Ph 


47.0150 


0.4335 


0.2050 


in 


11404 


Ph 


26.9260 


0.8671 


0.1950 




11697 


Ph 


0.1350 


-0.5553 


r\ n a on 

0.2480 




12110 


Ph 


338.6770 


0.0114 


0.2820 




12310 


Ph 


351.5230 


1.0208 


0.1510 




12326 


Ph 


326.2440 


1.0698 


0.2930 




12779 


Sp 


309.4730 


1.2194 


0.0800 


f 


12780 


Sp 


322.1570 


1.2302 


0.0490 


in 


12781 


Sp 


5.4080 


-1.0106 


0.0840 


f 


12825 


Ph 


359.5440 


-1.1823 


0.1500 


f 


12843 


Sp 


323.8790 


-0.9796 


0.1670 


f 


12856 


Sp 


332.8650 


0.7556 


0.1720 


in 


12897 


Sp 


18.4240 


-0.1027 


0.0170 




12903 


Ph 


327.6490 


0.3851 


0.0530 


HI 


12909 


Ph 


50.7700 


0.2366 


0.1310 




12917 


Ph 


17.0230 


1.2028 


0.2030 


111 


12923 


Ph 


9.9910 


0.6400 


0.1130 


f 


12930 


Sp 


309.6830 


-0.4764 


0.1470 


1 


12950 


Sp 


351.6670 


-0.8406 


0.0830 


r 


12964 


I'll 


309.9510 


-0.0693 


-0.0000 


f 


12971 


Sp 


6.6480 


-0.3033 


0.2350 




12979 


Sp 


11.6010 


0.0024 


0.1160 




12983 


Sp 


16.4580 


0.1454 


0.2650 




13007 


I'll 


4.8410 


0.6031 


0.1050 


r 


13013 


Ph 


14.3390 


0.8482 


0.1070 




13035 


I'll 


337.6120 


0.9692 


0.1310 




13063 


I'll 


338.9400 


-1.1845 


0.1340 


in 


13068 


I'll 


345.3330 


-0.6541 


0.0750 


f 


13070 


Sp 


357.7850 


-0.7466 


0.1990 


HI 



"Supcrnovac classification, Sp=Spectroscopic, Ph=Photomctric 
(see Section [2{ 

^Sample selection, f=final host spectroscopy sample (sec 
Section l2.3l l. m=Master sample (see Section 13.4.21 ) 
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Table Al - continued 



SN ID 


Class. a 


R,a 

( dc g) 


Dec 

(deg) 


H,easniit 


Sampli 


13071 


Ph 


358.6180 


-0.7188 


0.1790 




13099 


Sp 


359.8190 


-1.2507 


0.2660 




13113 


Ph 


19.7740 


-1.2338 


0.1220 


in 


13135 


Sp 


4.1740 


-0.4252 


0.1050 


f 


13254 


Sp 


42.0590 


-0.3468 


0.1810 


m 


13354 


Sp 


27.5650 


-0.8867 


0.1580 


f 


13458 


Ph 


16.4620 


-0.2496 


0.3190 




13511 


Sp 


40.6110 


-0.7942 


0.2380 




13522 


I'll 


21.6000 


-0.1613 


0.1670 




13545 


Ph 


52.3430 


0.5963 


0.2140 




13601 


Ph 


16.4310 


-0.5327 


f\ O A A f\ 

0.2 110 




13610 


Sp 


326.0140 


0.7255 


0.2980 




13633 


Ph 


4.6650 


0.0055 


0.3880 




13998 


Ph 


324.0760 


-0.0256 


0.1230 




14137 


Ph 


56.5190 


i \ a r\f\c\ 

-0.4009 


0.3120 




14153 


Ph 


49.3390 


0.0785 


0.1820 


f 


14193 


Ph 


359.0010 


1.0131 


0.1510 




14269 


Ph 


54.7590 


0.2792 


0.2810 




14279 


Sp 


18.4880 


0.3714 


0.0450 


f 


14284 


Sp 


49.0490 


-0.6010 


0.1810 


in 


14318 


Sp 


340.4250 


-0.1369 


0.0580 


f 


14340 


Ph 


345.8270 


-0.8553 


0.2780 




14377 


Sp 


48.2640 


-0.4718 


0.1390 


in 


14398 


Pll 


50.7100 


0.2729 


0.1180 




14421 


Sp 


31.8300 


1.2520 


0.1750 


f 


14753 


I'll 


44.3410 


0.6866 


0.1350 




14816 


Sp 


336.7160 


0.5064 


0.1070 


f 


14844 


Pll 


353.4150 


0.1766 


0.4750 




14961 


I'll 


15.9190 


0.9313 


0.3710 




15129 


Sp 


318.9020 


-0.3217 


0.1980 


f 


15136 


Sp 


351.1620 


-0.7179 


0.1490 


f 


15161 


Sp 


OF" 1 > \ ' \ ( \ 

35.8430 


0.8190 


0.2500 




15222 


Sp 


2.8520 


0.7020 


0.1990 


m 


15234 


Sp 


16.9580 


r\ o o o n 

0.8286 


0.1360 


in 


15262 


Pll 


342.4330 


-0.4042 


0.2470 




15421 


Sp 


33.7410 


0.6027 


0.1850 


m 


15425 


Sp 


55.5610 


0.4784 


0.1600 


f 


15443 


Sp 


49.8670 


-0.3180 


0.1820 


in 


15454 


Pll 


327.8590 


-0.8492 


0.3830 




15456 


Sp 


331.8630 


-0.9029 


0.0890 


in 


15467 


Sp 


320.0200 


-0.1774 


0.2100 




15533 


Pll 


353.8630 


-0.3733 


0.3440 




15568 


Pll 


356.8460 


1.0501 


0.2070 




15587 


Pll 


FT A A T OA 

54.4180 


f\ i it \ O A 

0.9984 


0.2190 




15648 


Sp 


313.7190 


-0.1958 


0.1750 


f 


15727 


Pll 


351.4000 


-0.0145 


0.1050 




15753 


I ) u 

rn 


/ .41 (\i 


-0.02 ( ( 


0.1590 




15765 


I'll 


32.8480 


0.2462 


0.3050 




15784 


Ph 


356.6760 


-0.6167 


0.2770 




15950 


Ph 


6.1890 


0.9841 


0.2200 


f 


15971 


Ph 


40.1130 


0.5263 


0.3160 




16038 


Ph 


55.5710 


-0.2378 


0.4200 




16069 


Sp 


341.2450 


-1.0064 


0.1290 


f 


16099 


Sp 


26.4210 


-1.0546 


0.1970 


in 


16163 


Pll 


31.4990 


-0.8558 


0.1550 


f 


16172 


Pll 


57.7050 


-0.2182 


0.2190 




16211 


Sp 


348.1630 


0.2660 


0.3110 




16215 


Sp 


18.4070 


0.4237 


0.0470 





"Supernovae classification, Sp=Spectroscopic, Ph=Photomctric 
(see Section [2]) 

''Sample selection, f=final host spectroscopy sample (sec 
Section l2.3l l. m=Master sample (see Section 13.4.21 ) 
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Table Al - continued 



SN ID 


Class. a 


R,a 

( dc g) 


Dec 

(dcg) 


H.easniit 


Sampl 


16287 


Sp 


46.6650 


0.0620 


0.1070 


m 


16259 


Sp 


352.0330 


0.8582 


0.1190 


111 


16276 


Sp 


nn t / • / 1 

20.5760 


1.0082 


0.0560 


111 


16280 


Sp 


1/11 o nn 

14.122U 


-1.2268 


0.0380 


I 


16314 


Sp 


onn t \<~m \ t\ 

320.9290 


-0.8430 


n / 1 / 1 • > n 

0.0630 


111 


16333 


Sp 


ooo nn a n 

328.9940 


-1.0703 


n nTnn 

0.0720 


111 


16392 


Sp 


ot n p" on 

27.9520 


0.2638 


0.0590 


I 


16462 


Ph 


1 t n 1 1 n 

17.0410 


-0.3859 


0.2450 




16482 


Sp 


o o o Tn o n 

328.7080 


n n n t 

0.9307 


0.2110 




16543 


Ph 


ft a on /in 

34.3940 


0.6481 


n nn 

0.2890 




16666 


Ph 


F o n o nn 

52.0300 


0.1031 


0.1330 


111 


16692 


Sp 


o on o o o n 

320.3830 


n nn /i n 

0.9949 


0.0340 


1 


16789 


Sp 


/i o to /i n 

43.7840 


0.2335 


0.3250 




17117 


Sp 


/in n on 

40.6020 


-0.7952 


0.1400 


I 


17134 


Sp 


54.5680 


-0.1101 


n no Tn 

0.0870 




17135 


Sp 


56.5290 


0.3903 


0.0310 




17171 


Sp 


326.5030 


-1.2193 


0.1600 


I 


17176 


Sp 


334.4030 


0.6133 


0.0930 


111 


17186 


Sp 


31.6160 


-0.8981 


0.0800 


111 


17206 


I'll 


45.9860 


0.7282 


0.1560 


111 


17215 


Sp 


54.9270 


1.0927 


0.1810 


111 


17219 


I'll 


3.1940 


-0.1031 


0.1940 


111 


17240 


Sp 


8.6420 


-1.2160 


0.0730 


I 


17280 


Sp 


55.7920 


0.1040 


0.1310 


I 


17316 


Ph 


3.7660 


-0.6243 


0.1040 




17332 


Sp 


43.7720 


-0.1477 


0.1830 


111 


17340 


Sp 


41.2130 


0.3653 


0.2570 




17366 


Sp 


o i r" to p" n 

315.7850 


-1.0312 


0.1390 


111 


17376 


Ph 


330.6440 


0.6728 


0.0720 




17434 


Ph 


1 o zioon 

18.4380 


n n to n 

-0.0730 


0.1790 


111 


17497 


Sp 


37.1360 


-1.0428 


0.1450 


111 


17458 


Ph 


o p" o i p" n n 

358.1590 


n n 

-0.3380 


0.0810 




17473 


Ph 


12.2590 


0.5474 


0.1700 


111 


17547 


Ph 


o o n i t • i n 

339.1610 


1.0913 


0.0610 




17605 


Sp 


309.2030 


0.0985 


0.1460 


I 


17629 


Sp 


on /•')/• n 

30.6360 


-1.0899 


0.1370 


I 


17784 


Sp 


52.4620 


0.0545 


0.0370 


111 


17880 


Sp 


/i /i nT /i n 

44.9740 


1.1601 


n n t n 

0.0730 


111 


17886 


Sp 


54.0070 


1.1048 


0.0410 


I 


17907 


Ph 


p* / > on Pf n 

56.2950 


n a nn a 

0.4094 


0.1940 


111 


17958 


Ph 


34.4350 


-0.7129 


0.2760 




18030 


Sp 


4.9330 


-0.4001 


0.1560 


111 


18047 


Ph 


oo nT /i n 

22.0740 


-0.6586 


0.3590 




18100 


Ph 


on a o on 

29.4320 


n A f\ A 

0.4243 


n a n n 

0.4980 




18201 


Ph 


47.3120 


-0.6450 


0.2930 




18224 


Ph 


O A O -1 rrf A 

348.1750 


-0.3129 


n nn 

0.3390 




18273 


I ) u 

rn 


334. (390 


-0.6303 


0.3160 




18298 


Sp 


18.2670 


-0.5400 


0.1200 


I 


18454 


Ph 


57.7460 


-0.2836 


0.3480 




18612 


Sp 


12.2880 


0.5966 


0.1150 


I 


18630 


I'll 


347.9800 


-0.2641 


0.3590 




18697 


Sp 


11.2240 


-0.9969 


0.1070 


I 


18721 


Sp 


3.0780 


-0.0777 


0.4030 




18751 


Sp 


5.7220 


0.7759 


0.0710 


I 


18764 


Ph 


4.0140 


1.1744 


0.2520 




18801 


Ph 


49.7600 


-0.4582 


0.2450 




18809 


Sp 


50.8810 


0.6673 


0.1320 


I 


18835 


Sp 


53.6850 


0.3555 


0.1230 


f 



"Supcrnovac classification, Sp=Spectroscopic, Ph=Photomctric 
(see Section [2{ 

b Samplc selection, f=final host spectroscopy sample (sec 
Section l2.3l l. m=Master sample (see Section 13,4.21 1 



Table Al - continued 



SN ID 


Class. a 


R,a 

(dcg) 


Dec 

(deg) 


H,easniit 


Sampli 


18855 


Sp 


48.6340 


0.2689 


0.1280 


m 


18890 


Sp 


1 r< a a on 

16.4430 


-0.7595 


n n/ -1 n n 

0.0660 


I 


18903 


Sp 


12.2510 


-0.3233 


0.1560 


I 


18927 


Sp 


A C! c a n 

46.684U 


-0.7573 


0.3610 




18959 


Sp 


r* a n on 

36.4080 


n Tn on 

0.7089 


0.4010 




19090 


I'll 


O P" T 1 c on 

357.1620 


-0.4065 


0.3120 




19155 


Sp 


01 n p" n 

31.2650 


0.1751 


n n TTn 

0.0770 


I 


19317 


I'll 


1 n ^ r rn 

310.4550 


1.0649 


0.1790 




19353 


Sp 


a 11 on 

43.1130 


0.2517 


0.1540 


111 


19616 


Sp 


t 1 nnn 

37.1000 


0.1860 


0.1650 


I 


19626 


Sp 


or noon 

35.9280 


-0.8265 


0.1130 


I 


19681 


I'll 


on a 1 on 

20.4120 


n a Tno 

-0.4702 


0.3510 




19778 


I'll 


a n r 1 n 

349.6880 


-0.4976 


n c\c n 

0.3960 




19787 


I'll 


0.2810 


n nnTn 

-0.0979 


0.1970 


I 


19794 


Sp 


f n 1 nn 

359.3190 


0.2485 


n on Tn 

0.2970 




19968 


Sp 


24.3490 


-0.3117 


0.0560 


I 


19969 


Sp 


31.9100 


-0.3240 


0.1750 


I 


20047 


I'll 


326.5910 


0.6311 


0.3740 




20064 


Sp 


358.5860 


-0.9172 


0.1050 


111 


20084 


Sp 


347.9770 


-0.5791 


0.0910 


I 


20141 


I'll 


357.5410 


-0.5245 


0.3410 




20314 


I'll 


4.1790 


0.7225 


0.2100 




20331 


I'll 


7.5400 


1.2459 


0.1840 


111 


20350 


Sp 


312.8070 


-0.9578 


0.1290 


I 


20386 


I'll 


55.7980 


-0.4896 


0.3080 




20420 


Sp 


338.8710 


0.4822 


0.1510 




20480 


I'll 


357.2740 


0.9182 


0.1680 




20528 


Sp 


a 1 1 n 

43.1210 


-1.1394 


0.1360 


111 


20625 


Sp 


5.6830 


-0.4794 


0.1080 


I 


20626 


I'll 


8.4760 


-0.5930 


n ot/ - * n 

0.2760 




on/ -1 to 

20678 


I'll 


17.5050 


1.2353 


0.2060 




20721 


I'll 


000 1 p" n 

323.1850 


n r 1 00 t 

-0.6227 


0.2120 




20726 


I'll 


42.2800 


-0.1595 


0.3200 




OnTO T 

20787 


I'll 


50.5110 


-0.4425 


n Tnn 

0.2700 




20788 


I'll 


51.6730 


-0.4778 


0.3930 




ooq on 
20889 


Sp 


52.3810 


0.5168 


0.20y0 


111 


20979 


I'll 


44.7730 


-0.2587 


0.1270 




21034 


Sp 


28.1420 


1.2441 


0.1090 


I 


21062 


Sp 


333.4290 


0.3966 


0.1390 


111 


21442 


I'll 


6.9600 


0.6003 


0.2130 




21502 


Sp 


353.6000 


-0.8902 


0.0890 


I 


21510 


Sp 


7.3520 


0.8311 


0.1500 




21669 


Sp 


11.6140 


-1.0609 


0.1240 


I 


21709 


I'll 


326.5280 


-1.0493 


0.1590 


111 


21872 


I'll 


7.6130 


-0.5915 


0.2290 




22006 


I'll 


48.7830 


-0.9626 


0.3960 




22075 


Sp 


29.9640 


1.2166 


0.1300 


f 



"Supernovac classification, Sp=Spectroscopic, Ph=Photornctric 
(see Section [2{ 

6 Samplc selection, f=final host spectroscopy sample (see 
Section l2,3| l. m=Master sample (see Section 13. 4.211 
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